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1.0  INTRODUCTION  AND  BACKGROUND 

Our  long-range  goal  is  to  understand  the  micromechanics  of  fatigue 
crack  initiation  and  to  use  this  insight  to  construct  a  model  of  the  stochas¬ 
tic  variability  of  crack  initiation  under  the  spectrum  loads  pertinent  to 
aircraft.  We  have  focused  our  attention  on  situations  where  cracking  is 
caused  by  peak  loads  below  the  alloy  elastic  limit.  The  genesis  of  the  crack¬ 
ing  is  highly  localized  plastic  deformation,  usually  in  individual  gr3irs  of 
an  alloy  and  at  its  surface.  Exciting  and  important  discoveries  a*ojt  tne 
mechanical  properties  of  surfaces  have  been  made  by  this  research.  >pse  are 
essential  ingredients  in  the  mathematical  description  of  the  mi ; r-o-ne; "a'- cs  o* 
localized  deformation,  an  area  in  which  we  have  also  made  advances.  T^is 
report  documents  the  progress  which  has  been  made.  Tne  major  effort  has  been 
to  develop  theoretical  tools  and  models  needed  to  extract  the  mechanical  prop¬ 
erties  of  the  surfaces  from  experimentally  measured  local  microplastic 
strains.  These  same  models  are  the  core  of  our  representation  of  the  microme¬ 
chanical  behavior  of  the  surface  and  will  be  a  starting  place  for  any  future 
formulation  of  a  description  of  the  overall  stochastic  behavior  of  crack 
initiation  in  the  presence  of  surface  deformation.  A  parallel  program  dealing 
with  surface  microplasticity  has  been  supported  by  NSF.*  Where  pertinent, 
results  from  the  NSF  contract  are  also  reported  here  with  appropriate  acknowl¬ 
edgement.  As  the  NADC/NAVAIR  research  began,  we  relied  upon  the  NSF  experi¬ 
mental  data  to  test  the  NAVAIR/NADC  surface  deformation  models.  Then,  as  this 
program  advanced  to  spectrum  loads,  additional  experimental  data  were  obtained 
with  the  NAVAIR  funding.  The  NSF  research  now  involves  4%  Cu-Al  alloys. 

The  major  findings  of  our  research  pertain  to  the  alloy  A1  2219-T851, 
the  model  material  used  in  this  study.  However,  there  is  reason  to  believe 
that  our  conclusions  are  generic  to  all  high-strength  structural  aluminums, 
and  have  important  implications  to  other  structural  materials,  especially  to 
certain  steels.  We  have  discovered  that  the  flow  stresses  within  the  larger 
grains  of  2219  decrease  progressi vely  with  fatigue,  ultimately  falling  to  less 
than  30*  of  their  prefatigue  value.  Usually,  the  surrounding  surface  remains 
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elastic.  These  changes  occur  slowly,  and  over  short  time  intervals  the  defor 
mation  behaves  as  if  the  local  flow  stress  was  a  constant.  The  effect  of 
loading  sequence  on  the  mean  stress  within  a  microplastic  grain  is  simply 
determined  by  the  reaction  stresses  created  when  the  plastic  deformation  is 
confined  by  the  nearly  elastic  surrounding  matrix.  Herein  is  the  apparent 
source  of  many  load  sequence  effects  on  fatigue  lifetime. 

We  begin  our  report  by  describing  in  detail  the  micromechanical 
models  of  the  deformation  of  a  surface  which  have  been  developed  (Sec.  2). 
Their  application  to  analysis  of  the  mechanical  properties  of  the  surface  of 
A1  2219-T851  is  then  described,  and  the  implication  of  the  models  to  the  re¬ 
sponse  of  the  alloy  to  selected  load  sequences  is  discussed  (Sec.  4). 

Finally,  we  briefly  discuss  our  findings  in  light  of  the  opportunities  for 
future  research  needed  to  develop  an  impl ementable  initiation  lifetime  pre¬ 
diction  methodology  (Sects.  5  and  6). 
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2.0  DEFORMATION  MODELS 

Our  theory  and  analysis  of  localized  surface  deformation  comprises 
several  elements.  These  are: 

(i)  The  description  of  the  constraints  imposed  upon  a  deformed 
grain  by  the  surrounding  elastic  matrix; 

(ii)  treatment  of  the  observed  nonuni formity  of  the  strain  fields 
within  a  grain,  especially  the  tendency  of  the  center  of  the 
grain  to  become  softer  with  fatigue  relative  to  the  boundary 
regions; 

(iii)  the  modelling  of  observed  hysteretic  properties  of  the  total 
strain  as  a  function  of  the  applied  stress;  and 

(iv)  the  calculation  of  plastic  work  accrued  in  a  grain  during 
fatigue,  for  future  correlation  with  the  statistics  of 
particle  fracture  as  a  test  of  the  crack  initiation  modelling 
conjectures . 

These  theoretical  considerations  raise  very  complex  questions,  and 
involve  many  areas  in  which  the  available  information  is  insufficient  to  de¬ 
cide  with  certainty  what  is  the  best  or  even  a  valid  model.  On  many  points, 
the  experimental  results  reported  in  Sec.  4.0  are  the  only  quantitative  data 
available,  representing  the  limit  of  guidance  to  theoretical  modelling. 
Therefore,  we  have  chosen  to  use  the  simplest  approach  possible  at  every 
stage,  not  diverting  effort  to  details  and  exactitudes  that  are,  in  our  best 
judgement,  most  likely  ancillary  to  the  determining  physics.  We  have  adopted 
a  phenomenological  approach,  conjecturing  a  reasonable  form  for  a  given  func¬ 
tion,  and  then  determining  various  parameters  by  direct  comparison  with  the 
data.  By  this  approach  we  have  allowed  the  experimental  results  to  speak  for 
themselves  as  far  as  possible,  rather  than  imposing  a  priori  theoretical 
prejudices . 
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We  shall  begin  in  Sec.  2.1  by  describing  the  theory  of  homogeneous 
deformation  in  a  single  grain.  Expressions  are  derived  for  the  plastic  and 
total  strains  as  a  function  of  the  applied  load.  These  are  used  in  Sec.  4.0 
to  deduce  the  flow  stress  and  strain  hardening  parameter  (or  plastic  stiffness 
constant)  from  experimental  hysteresis  loops.  Various  models  for  the  yield 
stress  and  strain  hardening  are  investigated.  The  models  for  the  yield  stress 
generally  correspond  to  von  Mises'^  criterion,  appearing  in  different  forms 
depending  on  whether  the  local  stress  state  is  assumed  to  be  uniaxial  or  bi¬ 
axial.  Strain  hardening  is  treated  either  by  Prager's  kinematic  hardening 
model,  or  by  a  simplified  hybrid  model  that  allows  certain  memory  effects  in 
the  local  stress-strain  response.  We  shall  also  describe  application  of  these 
models  to  loading  sequences  consisting  of  one  block  at  constant  amplitude  fol¬ 
lowed  by  another  at  reduced  amplitude.  It  will  be  shown  that  experimentally 
observed  memory  effects  manifest  in  the  transient  plastic  response  following 
the  change  in  load  amplitude,  can  help  to  determine  the  correct  model  for 
strain  hardening. 

A  model  of  certain  observed  inhomogeneities  in  the  strains  in  a  grain 
will  be  introduced  in  Sec.  2.2.  This  model  has  been  used  (Sec.  4)  to  analyze 
surface  strain  data  for  A1  2219-T851 ,  where  the  center  and  boundary  regions  of 
a  typical,  microplastic  grain  were  found  to  have  different  plastic  strains, 
and,  by  inference,  different  mechanical  properties.  This  inhomogeneity  is 
modeled  by  considering  the  grain  boundary  region  to  be  a  subdomain,  homoge¬ 
neous  within  itself,  nested  in  the  rest  of  the  grain  which  is  also  homogene¬ 
ous.  This  "dual  domain"  model  is  able  to  simulate  an  extraordinary,  and 
experimentally  observed,  anticlockwise  motion  in  the  hysteresis  loops  of  the 
measured  total  strain  vs  external  stress.  This  phenomenon  implies  that  the 
center  of  the  grain  suffers  extreme  softening  during  fatigue. 

In  Sec.  2.3,  a  generalization  of  the  dual  domain  model  is  presented 
that  deals  with  a  biaxial  stress  state,  which  is  the  general  case  for  a  grain 
at  a  free  surface.  We  believe  that  this  model  will  be  a  powerful  tool  for 
studying  not  only  the  biaxial  stress  state,  but  also  all  kinds  of  hysteresis 


w-w  *rw  w\;  »fYW  < 


r*^f  v^^’.v,ry^yT \r»y 


NADC-87046-60 


effects.  A  computer  code  implementing  the  model  has  been  completed.  It  is 
the  means  to  an  eventual  understanding  of  damage  accumulation,  and  the  proto¬ 
type  of  a  functioning  deformation  model  valid  under  all  kinds  of  spectrum 
loading . 

For  various  reasons,  including  the  general  state  of  theoretical  opin¬ 
ion  and  the  import  of  certain  previously  published  results  on  particle  frac¬ 
ture,  the  total  plastic  work  done  may  be  a  good  measure  of  the  status  of 
fatigue  damage  leading  to  initiation.  Therefore,  in  Sec.  2.4  we  discuss 
briefly  the  accumulation  of  plastic  work  in  a  homogeneous  grain.  We  show  how 
this  would  be  readily  included  in  the  computer  program  for  the  dual  domain 
model,  and  illustrate  the  importance  of  further  careful  studies  of  hysteresis 
and  memory  effects  in  deciding  the  aptness  of  assumed  work  hardening 
behaviors . 

At  the  end  of  this  section,  we  present  a  brief  review  (Sec.  2.5)  of 
other  various  modelling  approaches  we  have  considered  over  the  course  of  this 
program.  Although  none  of  these  approaches  is  per  se  part  of  the  set  of 
models  described  in  Secs.  2. 1-2. 4,  they  nevertheless  provided  decisive  and 
valuable  information  at  various  stages  of  our  learning  process.  Several  of 
them  remain  outstanding  candidates  for  more  detailed  theories  in  the  future. 

2 .1  The  Deformation  of  a  Homogeneous  Grain  Under  Uniaxial  Stress 

There  are  important  cases  where  the  inhomogeneity  of  a  grain  must  be 
taken  into  account,  as  will  be  described  later.  However,  in  many  experimental 
observations  the  grains  apparently  deform  uniformly,  and  even  in  an  inhomogen¬ 
eous  grain,  most  of  its  interior  region  may  still  be  regarded  as  homogeneous. 
To  analyze  hysteresis  data  for  all  such  cases,  relationships  between  the  plas¬ 
tic  strain  and  the  local  and  external  stress  were  derived  for  a  homogeneous 
grai n . 

The  grain  was  assumed  to  be  ellipsoidal  and  embedded  in  a  wholly 
elastic  matrix.  The  effect  of  the  free  surface  (the  specimen's  surface)  was 
neglected  because  of  numerical  demonstrations  of  its  relative  unimportance,  as 
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described  in  Sec.  2.5.  The  material  was  assumed  to  be  elastically  isotropic, 
and  Young's  modulus  Ee  was  taken  to  be  uniform  everywhere. 

According  to  Eshelby,^  if  the  plastic  deformation  in  an  ellipsoidal 
inclusion  is  homogeneous,  then  the  internally  constrained  strain  and  stress 
are  also  homogeneous.  In  this  sense,  the  assumption  of  homogeneity  in  the 
plastic  strain  is  self-consistent,  and  one  may  easily  write  down  the  local 
stress  in  terms  of  the  applied  load  and  the  state  of  plastic  strain. 

There  is  no  evidence  of  strong  slip  bands  in  the  measured  strains  nor 
of  shearing  deformation  of  the  overall  shape  of  the  observed  grains.  There¬ 
fore,  it  is  assumed  that  the  plastic  deformations  are  pure  extensions  or  con¬ 
tractions  in  the  frame  of  reference  defined  by  the  axes  of  the  ellipsoid.  In 
other  words,  the  ellipsoidal  axes  coincide  with  the  principal  axes  of  the 
plastic  strains.  As  illustrated  by  Fig.  2.1,  we  will  consider  one  axis  of 
symmetry  the  ellipsoidal  grain  to  be  parallel  to  the  x-axis  which  is  also  the 
stress  axis  for  the  uniaxial  external  stress  state.  If  we  then  write  for  the 
components  of  the  plastic  strain  cp  that 


we  may  consider  the  state  of  plastic  strain  to  be  entirely  defined  by  exp  and 
the  ratios  p  and  q.  The  plastic  strains  are  assumed  to  conserve  volume,  so 
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Fig.  2.1  Schematic  illustration  of  a  cross-section  of  the  ellipsoidal  grain, 
showing  strain  components.  The  x-axis  is  parallel  to  the  external 
stress  and  z  is  taken  to  be  normal  to  the  surface. 


The  components  of  the  local  stress  ox  may  be  written  in  terms  of  the  applied 
stress  oa  and  tp.  In  particular, 


where  a  is  a  function  of  p  and  q  and  depends  on  the  shape  of  the  grain.  It  is 
derived  from  Eshelby's  tensor  Sjjk-|,  which  relates  the  total  constrained 
strain  to  the  unconstrained  transformation  strain  in  an  ellipsoidal  inclusion. 
Parameter  a  usually  lies  between  0.5  and  1  (see  Sec.  2.2  for  further  discus¬ 
sion  and  Fig  .  2.7) . 


2.1.1  Stress  State 

For  the  simplest  models  in  which  the  stress  state  is  always  assumed 
to  be  approximately  uniaxial,  Eq.  (3)  is  the  only  statement  needed  to  estab¬ 
lish  the  relationship  between  ax,  aa  and  ep.  If  the  local  stress  is  to  be 
treated  as  biaxial,  then  the  other  components  of  ax  and  ep  must  be  consid- 
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ered.  Discussion  of  the  biaxial  case  is  deferred  to  Sec.  2.3.  We  shall  first 
consider  in  detail  the  behavior  of  models  based  on  the  assumption  that  the 
stress  state  is  uniaxial.  This  turns  out  to  be  approximately  correct  even  if 
the  biaxial  stress  is  included,  and  permits  us  to  derive  some  very  enlight¬ 
ening  expressions  for  certain  important  characteristics  of  the  observed 
hysteresis  loops. 

Even  on  the  macroscopic  scale,  where  experiments  are  much  easier, 
there  remains  doubt  and  controversy  over  the  best  way  to  describe  hysteresis 
effects  arising  from  plastic  flow  under  various  load  spectra.  On  the  micro¬ 
scopic  scale  of  the  interior  of  a  grain,  no  one  knows  the  true  state  of 
affairs.  We  are  presenting  the  first  data  on  the  subject.  There  is  no  cur¬ 
rent  knowledge  of  the  manner  of  strain  hardening  under  monotonic  or  cyclic 
loads,  the  degree  or  existence  of  Bauschi nger* s  effect,  other  memory  effects 
in  the  shape  of  the  yield  locus,  or  the  effects  of  fatigue.  Therefore,  our 
approach  has  been  to  apply  simple  models  containing  just  three  parameters  to 
be  determined  by  fitting  to  experimental  data.  These  may  be  identified  with  a 
microplastic  flow  stress  oQ,  a  linear  strain  hardening  parameter  (or  plastic 
stiffness)  Ep,  and  a,  which  determines  the  reaction  stress  and  depends  upon 
grain  shape.  Limits  on  the  range  of  a  can  be  defined  from  the  surface  grain 
shape,  but  a  also  depend  upon  grain  depth.  However,  caution  must  be  used  in 
assessing  the  physical  significance  of  these  parameters  because  their  values 
are  to  a  degree  model -dependent .  We  will  discuss  this  further  later. 

It  is  assumed,  in  the  the  uniaxial  models,  that  the  local  stress  and 
plastic  strain  are  always  linear  functions  of  the  applied  stress.  In  the 
absence  of  incremental  plastic  strain,  this  is  clearly  valid.  During  plastic 
strain,  it  amounts  to  saying  that  the  strain-hardening  is  linear,  which 
results  in  the  equilibrium  hysteresis  loop  for  the  grain  always  being  a 
parallelogram.  By  all  our  knowledge  of  macroscopic  plasticity,  this  is  a 
reasonable  model  capable  of  describing  the  experimentally  observed  loops. 
However,  as  we  shall  discuss  in  Sec.  2.5,  the  implications  for  studying  the 
accumulation  of  plastic  work  during  fatigue  require  careful  thought. 
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2.1.2  Yield  Conditions 

Saint  Venant^  guessed  the  now  established  fact  (generalized  by  Levy**) 
that  increments  in  the  components  of  plastic  strain  are  proportional  to  the 
corresponding  components  of  the  deviatoric  stress.  In  the  deforming  grain,  we 
must,  of  course,  consider  the  local  stress.  If  this  is  uniaxial,  then  it  fol¬ 
lows  that  p  and  q  of  Eqs.  (1)  and  (2)  are  always  constant  and  equal  to  -1/2. 
Therefore,  the  factor  a  in  Eq.  (3)  is  also  constant,  and  it  also  follows  that 
the  linear  strain  hardening  conjecture  may  be  expressed  as 

(4) 

expression  for  the  increment  in  total  strain  in  the 


leading  to  the  following 
plastic  regime: 


de 


t 

x 


If  the  material  is  not  yielding,  then 


(5) 


dc 


t 

x 


(6) 


Equations  (3)-(6)  are  immediately  sufficient  to  determine  the  increments  in 
plastic  and  total  strain  given  an  increment  in  the  applied  load.  To  determine 
the  entire  hysteresis  loop,  a  yield  criterion  must  be  included  to  distinguish 
between  the  elastic  and  plastic  regimes. 


In  the  uniaxial  case,  the  yield  criterion  of  von  Mises,  that  the 
second  invariant  of  the  deviatoric  stress  tensor  reach  a  critical  value,  is 
simply 
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or 

oj  =  /I  k  ,  (7) 

where  primes  appearing  as  prefixes  denote  a  deviatoric  quantity.  The  quantity 
/J  k  is  the  yield  stress  under  uniaxial  tension,  and  we  denote  it  oQ  in  our 
models.  (However,  as  remarked  above,  aQ  is  strictly  a  parameter  of  the  model, 
and  must  be  interpreted  carefully.) 

The  presence  or  absence  of  Bauschi nger ' s  effect  is  anticipated  by 
considering  two  models  of  the  way  the  yield  point  in  one  loading  reversal  is 
changed  by  the  magnitude  of  the  maximum  plastic  strain  achieved  during  the 
prior  loading  reversal.  One  model,  which  we  shall  call  the  "stationary" 
model,  assumes  that  oQ  is  unchanged  by  plastic  strain.  This  is  clearly  valid 
only  for  fully  reversed  cyclical  loading,  but  this  is  the  case  in  all  the 
experiments  to  be  reported.  This  criterion  may  be  written 


s  i 

do  •  o  >  0 
x  x 

and 


(3) 


where  oQ  is  constant.  In  the  "stationary"  model,  the  only  memory  of  prior 
loading  (such  as  Bauschi nger' s  effect),  arises  from  the  fact  that  the  local 
stress  in  a  grain  is  affected  by  the  constraint  of  its  plastic  deformation. 
The  other  yielding  model,  which  we  shall  call  the  "kinematic"  model,  is 
Prager's^  model  of  kinematic  hardening.  Following  plastic  deformation  during 
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one  reversal,  flow  in  the  opposite  direction  in  the  next  reversal  always 
occurs  when  the  magnitude  of  the  local  stress  changes  by  2o0  from  its  final 
value  0-j  on  the  prior  reversal.  If  n  =  +  1  denotes  the  sense  of  the  current 
reversal ,  the  "kinematic"  criterion  for  yield  is 

n(  ax  *  °i )  >  2  °0  *  ( 9 ) 


2.1.3  Equi  1  i  b riu m  Loo p  Re sponse 

In  fully  reversed  loading,  the  properties  of  these  deformation  models 
verifiable  by  experiment  are  easily  discovered  with  some  simple  algebra.  The 
experiments  involve  the  measurement  of  total  local  strains  (c^)  or  loop  widths 
(W)  vs  the  applied  external  stress  0a.  The  following  results  neglect  the 
biaxial  stress  components  induced  by  the  localized  plastic  deformation.  Under 
fully  reversed  cyclic  loading  of  constant  amplitude,  both  yielding  models  pre¬ 
dict  that  the  total  strain  vs  external  load  hysteresis  will  be  a  parallelo¬ 
gram.  Since  this  should  correspond  to  the  observed  hysteresis  loop,  we  shall 
consider  it  in  detail  . 

Let  the  maximum  amplitude  of  the  applied  stress  be  s,  and  let  the 
onset  of  yield  during  the  loading  half  of  a  cycle  occur  when  0a  =  S}  (see  Fig. 
2.2).  Then  the  characteri sties  of  the  hysteresis  loop  are  as  follows: 

( a )  "Stationary  Model " 


S1  = 


20r)(EP  +  oEe)  -  oEPs 
2EP  +  aEe 


o0  +(o0  -  s) 


aE 


2EP  +  oE6 


S1  4  °o 


(10) 


11 


»  J 


» •  .1 


NADC-87046-60 


Fig.  2.3  The  width  of  the  hysteresis  loop  of  total  strain  vs  applied 

load  at  zero  applied  load  for  the  "stationary"  model,  (a)  vs  og 
for  constant  s  and  Ep;  (b)  vs  s  at  constant  oQ  and  E  . 

closes  again  as  the  local  flow  stress  falls  from  s  to  zero.  The  curve  is 
symmetric  about  its  peak,  and  passes  through  the  origin  if  the  material  is 
el astic-perfectly  plastic  (Ep  =  0).  The  reduction  of  WQ  for  small  oQ  occurs 
simply  because  yield  occurs  before  the  external  load  returns  to  zero  on  each 
reversal;  i.e.,  for  s^  <  0.  On  the  other  hand,  the  maximum  plastic  strain 
continues  to  increase  as  a  decreases  (Eq.  (11)).  Note  that,  since  the  maxi 
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mum  loop  width  at  fixed  load  also  varies  inversely  with  the  plastic  stiffness 
Ep,  an  observed  decrease  in  loop  width  may  be  attributed  either  to  a  falling 
yield  stress  o0  or  a  rising  plastic  stiffness  Ep.  Additional  analysis  to 
determine  both  Ep  and  oq  and  the  source  of  observed  loop  width  closing  is 
considered  in  Sec.  4.0.  Figure  2.3b  shows  that,  for  fixed  o0.  the  loop  width 
is  monotonically  increasing  with  increasing  cyclic  stress  amplitude.  In  the 
perfectly  plastic  case  (Ep  =  0),  the  loop  width  remains  constant  when  the  load 
exceeds  2o0. 

( b )  "Kinematic  Mooel " 

S1  =  2o0  -  s- 

p  =  s  ~  ao  4 

Sax  Ep  + 


(13) 

(14) 


t 

cmax 


(s 


1  -  a 

EP  +  oEe 


(15) 


The  width  WQ  of  the  loop  when  the  applied  load  is  zero  is  shown  as  a  function 
of  aj  in  Fig.  2.4a,  and  as  a  function  of  s  in  Fig.  2.4b.  WQ  in  the  "kine¬ 
matic"  model  does  not  decrease  with  falling  flow  stress,  but  remains  constant 
when  oq  <  0.5  s  (Fig.  2.4a).  For  constant  oq»  W0  is  a  strictly  monotonically 
increasing  function  of  s  (Fig.  2.4b).  Above  s  =  2cjq,  Wq  is  a  constant. 


2.1.4  Transient  Deformation  Behavior 

It  is  important  to  bear  in  mind  that  there  is  no  sound  a  priori  rea¬ 
son  to  prefer  either  the  "stationary"  or  "kinematic"  model.  In  fact,  we 
should  expect  that  the  truth  lies  somewhere  between  them,  i.e.  that  they  may 
perhaps  serve  as  limiting  cases  in  the  spectrum  of  plausible  hardening  models. 
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(a) 


0  °o  2  oo  s 

(b) 

Fig.  2.4  The  width  of  the  hysteresis  loop  of  applied  load  vs  total 

strain  at  zero  applied  load  for  the  "kinematic"  model,  (a)  vs 
eg  at  constant  E^.  (b)  vs  s  at  constant  E  . 


One  of  the  essential  differences  between  the  two  models  is  the  presence  or 
absence  of  memory  effects.  We  shall  now  demonstrate  this  by  considering  a 
load  sequence  comprising  a  first  block  of  fully  reversed  cycles  at  one  con¬ 
stant  amplitude  s  followed  by  another  at  a  new,  reduced  amplitude  s'  that  is 
still  above  the  flow  stress  og.  We  are  interested  in  the  existence  and  nature 
of  any  transient  behavior  at  the  beginning  of  the  second  loading  block. 


Consider  first  the  "kinematic"  model.  The  plastic  strain  at 
s'  >  on  >  Si  when  the  load  is  increasing  is 
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:P(s')  = 


-  er  + 
max 


EP  +  aEe 


S1  -  a 
_ 0 

EP  +  aEe 


(16) 


But,  comparison  with  Eq .  (14)  shows  that  this  is  just  the  maximum  strain  that 
would  be  experienced  in  equilibrium  if  the  load  amplitude  was  s'.  Therefore, 
when  the  load  amplitude  is  reduced  from  s  to  s' ,  the  material  at  once  follows 
the  equilibrium  loop  for  the  new  amplitude.  No  memory  of  the  amplitude  s  or 
the  existence  of  a  prior  higher  amplitude  block  is  present. 

This  is  not  the  case  for  the  "stationary"  model.  In  the  "stationary" 
model,  the  plastic  strain  at  s’  >  og  >  s^  is 


n  s  -  on  s  -  s 

£P  (  S  '  )  = - +  - 

X  2EP  ♦  aE6  EP  ♦  aEe 


E  (s-s’ ) 


2EP  +  aEe  ( 2EP  +  aEe)(EP  +  a£( 


(17) 


Comparison  with  the  maximum  plastic  strain  ^ax($')  at  equilibrium  with  a 

stress  amplitude  of  s'  (Eq.  (11))  shows  that  ep  (s')  is  less  than  ep  (s')  by 

x  roax 

the  second  term  of  Eq .  (17).  That  is  to  say,  at  any  point  s'  of  the  hyster¬ 
esis  loop,  the  deformation  suffered  on  the  prior  cycle  has  not  been  suffi¬ 
ciently  cancelled  and  reversed  to  be  at  the  equilibrium  level  ep  (s’)  corres- 

max 

ponding  to  s'.  (The  second  term  of  Eq.  (17)  is  always  greater  than  or  equal 
to  zero.  Equality  occurs  in  the  perfectly  plastic  limit  Ep  =  0.)  If  the  new 
block  at  amplitude  s'  begins  on  a  tensile  reversal  (see  Fig.  2.5),  then  we  may 
expect  the  following.  The  loop  width  measured  between  the  first  two  zeroes  of 
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TIME 


Fig.  2.5  Simulated  loading  history  for  the  study  of  transient  behavior. 

the  new  block  (points  A  and  B  in  Fig.  2.5)  will  be  greater  than  WQ(s')  for  the 
equilibrium  loop  at  load  s'.  If  s1  >  0,  it  is 


w:  =  w  (s')  +  - 


a(l  -  a) ( s  -s'  ) 


( E1^  +  aE  )  (2E*>oE  ) 


The  width  measured  over  the  compressive  half  cycle  between  the  next  two  zeroes 
(points  B  and  C)  can  be  shown  to  be  less  than  WQ(s').  The  loop  widths  for 
succeeding  reversals  approach  W0(s')  asymptotically,  and  numerical  examples 
show  that  the  new  equilibrium  loop  is  achieved  within  about  five  cycles. 

The  fact  that  the  "stationary"  model  shows  memory  in  this  transient 
region  while  the  "kinematic"  model  does  not  will  be  used  in  the  analysis  of 
data  (Sec.  4.0)  to  argue  against  Prager's  purely  kinematic  hardening  model. 
However,  it  should  not  then  be  inferred  that  the  "stationary"  model  is  there¬ 
fore  correct.  In  terms  of  the  total  strain  vs  local  stress  within  a  grain, 
the  "stationary"  model  shows  no  Bauschinger  effect,  which  may  well  lead  to 
irreconcilable  differences  with  other  experiments.  The  truth  probably  lies 
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between  these  two  extremes.  One  of  the  important  tasks  for  future  work  is  to 
devise  measurements  to  test  various  strain  hardening  conjectures,  under  both 
uniaxial  and  biaxial  stress  states.  This  is  a  prerequisite  to  the  development 
of  a  reliable  model  of  microplastic  deformation  during  stochastic  loading. 

2  .2  Anomalous  Hysteresis  in  Inhomogeneous  Grains  -  The  "Dual  Domain" 

Model 

For  the  sake  of  completeness,  the  following  description  of  the  dual 
domain  model  under  uniaxial  loading  is  presented.  This  model  was  formulated 
under  the  support  of  NSF  funding  when  it  became  clear  that  a  micromechanical 
approach  we  were  trying  (Sec.  2.5.3)  as  a  means  to  derive  quantitative  values 
for  oq  was  not  going  to  be  fruitful.  This  analysis  helps  one  visualize  the 
underlying  physics  of  the  deformation  which  leads  to  the  unusual  grain  bound¬ 
ary  stress-strain  response.  However,  the  actual  values  of  oQ  found  (Sec.  4.0) 
have  been  obtained  by  the  numerical  code  developed  for  NADC/NAVAIR  and  dis¬ 
cussed  in  Sec.  2.3,  because  the  strain  behavior  when  both  domains  are  plastic 
can  be  extremely  complex.  Later  it  was  realized  that  an  even  simpler  assump¬ 
tion  of  homogeneous  deformation  could  be  profitable  and  this  led  to  further 
analysis  (Sec.  2.1)  and  data  collection  (Sec.  4.0),  thus  permitting  calcu¬ 
lation  of  og,  a  and  from  experimental  data.  This  work  and  all  the  other 
theoretical  results  described  here  were  developed  under  the  NADC/NAVAIR 
contract . 

One  of  the  most  extraordi nary  experimental  phenomena  observed  in  A1 
2219-T851  is  the  existence  of  anticlockwise  hysteresis  in  measurements  of 
total  strain  as  a  function  of  external  stress.  That  means,  in  some  parts  of  a 
grain,  the  total  strain  on  tensile  loading  can  be  greater  than  that  at  the 
same  external  stress  during  unloading.  This  corresponds  to  an  anticlockwise 
motion  about  the  hysteresis  loop,  rather  than  the  normal  clockwise  motion. 

This  anticlockwise  motion  is  generally  observed  in  a  small  region  adjacent  to 
the  boundary  of  a  microplastic  grain,  while  the  rest  of  the  same  grain  exhib¬ 
its  the  normal  clockwise  hysteresis.  Continuing  with  the  assumption  that  the 
local  stress  state  is  approximately  uniaxial,  we  shall  devise  a  simple  model 
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that  can  reproduce  all  of  the  observed  hysteresis  phenomena.  We  then  describe 
in  Sec.  2.3  a  generalization  implemented  on  a  computer  to  allow  considerati on 
of  biaxial  stress  states  and  random  loads. 

In  the  "dual  domain"  model  the  inhomogeneity  is  introduced  in  the 
form  of  an  ellipsoidal  subdomain  of  the  grain.  As  it  turns  out,  with  further 
approximations ,  the  model's  properties  are  shown  to  be  independent  of  the 
precise  position  (but  not  the  orientation)  of  the  subdomain.  However,  we 
shall  generally  regard  it  as  being  at  the  grain  boundary.  Its  axes  are  taken 
to  lie  parallel  to  those  of  the  grain,  but  its  shape,  determined  by  the  semi¬ 
axes,  is  not  restricted. 

Suppose  that  part  of  the  grain  lying  outside  the  subdomain  (call  this 

p 

domain  1)  suffers  a  uniform  deformation  c..(l)  and  the  subdomain  (or  domain 

’  J  p 

2)  suffers  a  different  uniform  deformation  e..(2)  .  We  may  regard  this  state 

P  '  J 

as  the  superposition  of  the  deformation  eT  . ( 1 )  acting  over  the  entire  grain 

D  p  '  J 

and  the  deformation  e . ^ (2)  -  e-j(l)  acting  over  the  subdomain  (domain  2). 
Because  of  the  constraining  matrix,  reaction  stresses  accompany  each  of  these 
deformations .  As  for  the  homogeneous  grain,  we  require  simple  expressions  for 
these  reaction  stresses  in  order  to  calculate  the  local  stress. 

Because  of  Eshelby's  result,^  we  may,  as  for  the  homogeneous  grain, 
consistently  assume  that  as  cyclic  loading  progresses,  the  subdomain  (domain 
2)  will  remain  homogeneous.  In  the  grain  itself,  the  same  would  be  true  if 
not  for  the  reaction  field  arising  from  and  external  to  the  deforming  subdo¬ 
main.  The  presence  of  this  field  implies  nonuniformity  of  the  local  stress, 
and,  therefore,  of  the  plastic  and  total  strains  in  the  grain.  However,  if 
the  subdomain  is  a  small  part  of  the  entire  grain,  then  except  very  close  to 
subdomain  boundaries,  the  reaction  stress  to  subdomain  deformation  will  be 
very  small  relative  to  the  reaction  stress  created  by  the  grain's  own  defor¬ 
mation.  (We  have  reached  this  conclusion  from  many  numerical  calculations 
based  on  Eshelby's''  expressions  for  the  field  outside  an  ellipsoidal  inclu¬ 
sion.)  Therefore,  we  neglect  it.  This  amounts  to  assuming  that  the  plastic 
deformation  in  the  grain  does  not  depend  on  the  deformation  in  the  subdomain, 
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a  condition  that  happily  simplifies  the  model's  algebra.  The  local  stress  in 

p 

the  grain  remains  entirely  determined  by  e..(l),  as  it  was  for  the  homogeneous 

’  J 

grain.  It  is  then  consistent  to  assume  that  the  grain  (domain  1)  also  remains 
homogeneous  throughout  the  fatigue  process. 

In  an  expanded  notation  in  which  the  argument  1  refers  to  domain  1 
and  the  argument  2  to  domain  2,  we  now  have  that 

o*(l)  »  oa  -  a(l)Ee^(l)  (19) 

and  c*(2)  =  aa  -  a(l)EecP(l)  -  a(2)Ee(eP(2)  -  e£(l))  .  (20) 


The  coefficient  a(2)  is  analogous  to  a(l),  being  derived  from  Eshelby's  tensor 
as  described  in  Sec.  2.1.  In  the  present  section,  we  shall  see  that 
the  relative  values  of  a(l)  and  a(2),  which  depend  on  the  shapes  of  the  grain 
and  the  subdomain,  have  a  controlling  influence  over  the  nature  of  hysteresis 
in  the  subdomain. 

As  in  Sec.  2.1,  we  shall  assume  linear  strain  hardening,  so  that  an 
equation  of  the  forni  of  Eq.  (5)  can  be  written  for  both  domain  1  and  domain 
2.  Young's  modulus  is  assumed  uniform  everywhere,  but  the  plastic  stiffness 
Ep  may  be  different  in  the  two  domains,  so  that,  during  plastic  flow, 


and 


dej(l)  =  da*(l)^-  +  4 


EP  ( 1 ) 


dcj(2)  =  dc^( 2)  +  -p^— ; 

x  x  Ee  E  (2) 


(21) 


(22) 
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To  begin  study  of  the  "dual  domain"  model,  the  yield  criterion  of  the  "sta¬ 
tionary"  model  of  Sec.  2.1  will  be  assumed  and  flow  stresses  og( 1 )  and  og(2) 
will  be  attributed  to  the  two  domains.  Results  concerning  the  qualitative 
characteristics  of  hysteresis  are  unaffected  by  this  choice  of  model  for 
yielding. 

Let  us  examine  the  conditions  for  anticlockwise  hysteresis.  Anti¬ 
clockwise  hysteresis  occurs  if,  during  the  tensile  loading  part  of  the  cycle, 
the  slope  of  the  external  stress  vs  total  strain  curve  after  yielding  exceeds 
that  before  yielding  (see  Fig.  2.6).  In  either  domain,  the  slope  Se  of  the 
applied  stress-total  strain  curve  before  yielding  is  Ee.  In  domain  1,  the 
slope  SP(1)  after  yielding,  from  Eqs.  (19)-(22)  is 


SP(1) 


r  EP( 1 )  *  a(l)E 
1  EP(1)  +  Ee 


e 


(23) 


K»»  !)«)( 
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Fig.  2.6  Illustrating  how  anticlockwise  hysteresis  can  occur. 
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In  the  subdomain,  the  slope  SP(2)  after  yielding  is 
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semi-axis  less  than  a  few  percent  of  the  other  semi-axes).  This  is  not  be¬ 
lieved  to  be  the  case  in  the  experiments  reported  below.  Furthermore,  since  a 
can  barely  exceed  unity,  the  ratio  of  the  slopes  Se  to  SP(1)  can  never  equal 
that  implied  by  the  experimentally  observed  magnitudes  of  negative  loop 
widths . 

For  anticlockwise  hysteresis  to  occur  in  domain  2,  Sp(2)  must  exceed 
Se.  This  may  very  easily  happen.  For  example,  when  domain  2  is  always  elas¬ 
tic  (Ep(2)  -*•»),  as  it  sometimes  appears  nearly  to  be,  the  ratio  of  Sp(2)  to 
Se  is  just 


SP(2)  _  EP(1)  +  a(l)Ee  (25) 

Se  EP( 1)  +  a(2)Ee 

This  exceeds  unity  whenever  a(l)  exceeds  a(2),  and  Fig.  2.7  shows  that  this  is 

generally  the  case  when  domain  2  is  more  spherical  than  the  entire  grain. 

This  situation  is  perfectly  feasible.  It  is  evident  that  anticlockwise 
hysteresis  may  also  occur  readily  under  similar  conditions  concerning  the 
relative  shapes  of  the  domains,  if  EP(2)  is  finite.  For  any  given  a(l)  and 
a(2),  the  condition  may  be  stated  as 

EP(2)  ,  (1  -  a(2)IEP(l)  ♦  (1  -  a( 1 ) ) a( 2 )Ee  .  „ 

- adr-'jni - _>0-  (26) 

For  the  usual  case  that  0  <  a(l)  <  1  and  0  <  a(2)  <  1,  the  numerator  is  posi¬ 
tive,  and  we  see  that  a(l)  must  exceed  a(2)  for  a  critical  value  of  EP(2)  to 
exist.  Furthermore,  when  0  <  a(l),  a(2)  <  1,  the  right  hand  side  of  Eq .  (24) 
is  a  strictly  monotoni cal ly  increasing  function  of  EP(2).  Since  it  is  also 
less  than  Se  when  Ep(2)  =  EP(1),  it  follows  that  the  critical  value  of  EP{2), 

p 

if  it  exists,  is  always  greater  than  E  (1);  i  .e.,  the  subdomain  must  be  rela¬ 
tively  hard  for  anticlockwise  hysteresis  to  occur. 
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Note  that  no  reference  has  been  made  yet  to  the  flow  stresses  og(l) 
and  og(2),  except  the  implicit  assumption  that  the  local  stresses  have  ex¬ 
ceeded  them  at  some  point.  The  criterion  for  negative  loop  width  is  entirely 
determined  by  the  plastic  stiffness  parameters  and  the  shapes  of  the  domains. 
The  flow  stresses  og(l)  and  o0(2)  appear  in  the  expressions  for  the  magnitude 
of  the  loop  widths.  Furthermore,  in  the  event  that  both  domain  1  and  then 
domain  2  yield  on  each  reversal,  the  flow  stresses  determine  the  topology  of 
the  hysteresis  loop.  An  example  is  illustrated  schematically  in  Fig.  2.8. 

The  point  A  marks  the  onset  of  yield  in  domain  1  during  the  tensile  loading 
half-cycle.  Between  A  and  B,  domain  2  remains  elastic,  and  if  a(l)  >  a(2),  the 

slope  of  the  curve  increases  as  shown  in  the  figure.  B  marks  the  onset  of 

yielding  in  domain  2  and,  if  Ep{2)  is  less  than  the  threshold  given  in  Eq . 

(26),  the  slope  of  the  stress-strain  curve  will  become  less  than  it  was  when 
both  domains  were  elastic.  This  is  the  case  shown.  When  the  half-cycle  is 
completed  and  unloading  begins,  the  material  response  will  descend  on  a  line 
parallel  to  the  original  elastic  loading  line.  However,  depending  on  whether 
the  load  reversal  occurs  at  C  or  C',  the  descending  elastic  line  may  or  may 
not  intersect  the  loading  segment  AB.  Which  case  occurs  depends  on  the  rela- 


Fig.  2.8  Possible  topologies  for  the  hysteresis  loop  in  a  relatively  hard 

subdomain.  I 
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tive  magnitudes  of  og(l),  o0(2),  the  load  amplitude,  as  well  as  upon  the  plas¬ 
tic  stiffnesses  Ep(l)  and  Ep(2).  For  the  loop  ABC,  the  width  measured  at 
equal  external  stress  during  loading  and  unloading  will  appear  positive  (con¬ 
ventional)  at  high  stresses  and  negative  (anticlockwise)  at  low  stresses. 
Observations  of  such  patterns  have  indeed  been  made  and  are  reported  in  Sec. 
4.0. 

In  fact,  as  will  be  seen  in  Sec.  4.0  the  "dual  domain"  model  has  a 
very  high  degree  of  success  in  explaining  the  first  available  measurements  of 
strain  inhomogeneity.  Not  only  are  the  qualitative  observations  accounted  for 
but  even  the  values  of  the  model's  parameters  giving  best  fit  to  the  data  look 
very  plausible  and,  where  applicable,  are  in  good  agreement  with  the  values 
deduced  from  the  much  more  reliable  modelling  of  the  deformation  in  homoge¬ 
neous  grains.  Howeve^,  one  must  be  wary  of  being  seduced  by  this  success  into 
inferring  too  much.  Other  assumed  geometrical  distributions  of  inhomogeneity 
may  give  equally  appealing  results,  and  we  do  not  yet  have  enough  data  to  know 
whether  another  model  would  be  more  or  less  realistic.  Furthermore,  we  have 
not  yet  begun  to  investigate  how  the  various  observed  hysteresis  effects 
depend  on  the  nature  of  the  stress  state.  It  could  easily  be  that  an  analysis 
of  inhomogeneity  in  a  biaxial  stress  field  will  reveal  qual itati vely  different 
mechanisms  for  generating  anticlockwise  hysteresis.  In  the  next  subsection, 
we  shall  present  a  model,  already  implemented  on  a  computer,  to  tackle  just 
such  questions. 


2.3  The  "Dua’  Domain"  Model  in  the  Biaxial  Case 


The  "dual  domain"  model  has  been  generalized  to  consider  the  gene¬ 
rally  biaxial  stress  state  that  exists  in  a  grain  at  a  free  surface.  Since 
there  is  no  longe"  a  direct  proportionality  between  increments  in  the  compo¬ 
nents  of  plastic  strain  and  an  increment  in  the  applied  load,  it  is  no  longer 
possible  to  derive  simple  expressions  for  the  characteristics  of  hysteresis 
as  in  the  preceding  subsections.  Instead,  we  have  developed  a  numerical 
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approach,  for  which  a  computer  program  has  now  been  written.  Since  it  will  be 
the  basis  of  future  work  on  the  subject,  we  shall  describe  it  here  fully. 

In  the  numerical  model,  the  status  of  the  system  is  defined  at  any 
time  by  the  applied  stress  and  the  x  and  y  principal  components  of  plastic 
strain  in  each  domain  (the  free  surface  being  the  plane  z  =  0) .  Given  an 
increment  6oa  in  the  still  uniaxial  applied  stress,  the  change  in  local  stress 
state  is  given  by 

6cx(l)  =  cc3  -  a  6e^(l) 

X  '  x  X  x v 

6cx(l)  =  a  6e°( 1 ) 
y  y  x 

oax'2)  -  6c*  -  a;  czPx{2)  -  «;  6^(1) 

*«£«>  ■  <*$  * «;  ■ ' ?  -%»>  •  (27> 

It  is  to  be  remembered  that  coefficients  a,  a'  a"  a  ,  a'  and  a"  depend 

x  x  x  y  y  y 

not  only  on  the  grain's  shape,  but  also  on  the  components  of  the  deviatoric 
local  stress  tensor.  They  are  derived,  as  before,  from  Eshelby's  tensor 
Si jki  .  However,  since  the  ratios  of  the  components  of  the  deviatoric  stress 
tensor  are  always  changing  (in  contradistinction  to  the  uniaxial  case)  the 
coefficients  must  be  evaluated  afresh  at  every  increment  of  load. 

The  increments  in  plastic  strain  are  assumed  to  be  given  by  the 
Prandl-Reuss  equations:®’^ 


J  -  Jf  V 
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<(» 

6E>» 

6'x(2> 


3'°*<» 

.  do*(l) 

2EP(1) 

oX(l) 

3^(1> 

.  do*(l) 

2EP( 1 ) 

5*(1) 

3*  °x(2) 

.  deA(2) 

2  E  P  ( 2 ) 

c*(2) 

3‘c*(2) 

y 

.  dcX(2) 

2EP(2) 

3A(2) 

(28) 


where  a  preceding  prime  refers  to  a  deviatoric  quantity,  and 

dc"(i)  (2°x(i)  ’  °^i))6o^i)  +  -  <(i))6aj(i) 

FiT)  (<(i)  -  o^(i)  )2  4  (7(i)7  4  (o^Ci))2 


From  Eqs.  (27)  and  (23)  one  has  at  once  explicit  expressions  for  the  plastic 

P  P  d 

strain  increments  fee  (1)  and  bex(2)  i°  terms  of  6o  ,  which  may  be  regarded  as 
simultaneous  ordinary  differential  equations  of  the  first  order.  They  may  be 
solved  quite  adequately  by  simple,  low  order  techniques.  In  fact,  a  first 
order  method  has  been  used  satisfactorily. 

To  begin,  we  have  coded  Prager’s  kinematic  hardening  rule  for  the 

C  0 

translation  of  the  yield  locus  in  (a  (i),  o  ( i ) )  space  in  each  domain.  With 

x  y 

von  Mises  criterion,  the  yield  locus  for  domain  i  is  given  by 


i 


0,  ..  , •  X,  .. 

ex(i) )  0) 
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where  (0°(i),  o^(i))  defines  the  origin  of  the  yield  locus  at  any  time,  as 
x  y  o 

given  by  Prager's  pin  and  rigid  frame  model.  During  yielding,  if  the  local 

stress  changes  by  (6o£(i),  6o£(i))  in  (o*(i),  o£(i)l  space,  then,  according  to 

x  y  n  x  n  y 

Prager's  model,  the  translation  (6ox(i),  60  (i))  of  the  origin  of  the  yield 
locus  is  given  by 


6o°(i)  =  1/D  (oy£(-i)  -  2o^(i  )  )2  6oJ(i)  +  (o^(i )  -  2o*(i  )  )(o*(i )  -  2o£(  1 )  )6o£(i  ) 

6oy(i)  =  1/D  Cc^(i  )  -  2Sf(i))(Sf(i)  -  2cy£(i )  )6cf(i  )  +  (o£(i)  -  2oy£(i))2  6oy£(i) 

where  D  =  (Sy£(i)  -  2o£;i)}2  +  (o£:i)  -  2o£(i))2 

and  o£;i)  =  o£(i)  -  oj(i)  .  (30) 


When  using  Prager's  model,  the  origins  (o°(i),  o°(i))  should  be  regarded  as 

x  y 

independent  variables  constituting  part  of  the  status  of  the  system  at  any 
time . 

As  with  the  uniaxial  approximat i on ,  one  finds  that  the  kinematic 
model  just  defined  contains  no  memory.  If  block  loading  is  simulated,  the 
system  falls  into  its  new  equilibrium  hysteresis  loop  on  the  first  cycle  of  a 
new  block.  As  we  have  already  said,  this  contradicts  experimental  results 
presented  below.  Therefore,  a  necessary  improvement  of  the  computer  simula¬ 
tion  will  be  to  include  a  better  representation  of  strain  hardening  and 
Bauschi nger ' s  effect  than  Prager's  kinematic  model. 

In  its  present  form,  the  chief  benefit  of  the  computer  program  is 
that  it  allows  rapid  simulation  of  the  plastic  response  to  a  biaxial  local 
stress  state.  Furthermore,  there  is  no  restriction  whatever  on  the  uniaxial 
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applied  load.  This  may  follow  any  designated  spectrum,  since  the  load  se¬ 
quence  is  simply  defined  as  a  sequence  of  reversals  of  any  amplitude.  Because 
of  the  simplicity  of  the  model,  it  would  be  quite  feasible  to  follow  many 
thousands  of  reversals  in  detail.  Thus,  the  program  gives  us  a  powerful  tool 
to  study  the  role  of  accumulated  plastic  strain  or  total  plastic  work  done  in 
determining  initiation  rates  under  diverse  loading  histories.  A  phenomenolog¬ 
ical  and  testable  distinction  is  drawn  between  such  quantities  in  the  interior 
of  a  grain  and  near  its  boundaries. 

2 .4  Plastic  Work  Done 

c  P 

For  a  given  local  stress  c  and  a  given  plastic  strain  e  ,  the  total 
irrecoverable  plastic  work  done  is 


w 


p 


(31) 


Since  there  is  good  reason  to  expect  plastic  work  to  be  a  useful  measure  of 
fatigue  damage,  we  shall  consider  now  some  simple  expressions  for  the  work 
done  in  a  homogeneous  grain  under  constant  amplitude  loading,  using  the  "sta¬ 
tionary"  and  "kinematic"  uniaxial  models  described  above  (Sec.  2.1).  The 
purpose  of  this  exercise  is  to  get  a  feeling  for  the  effect  upon  the  plastic 
work  done  for  each  model  assumed  for  the  yield  criterion. 

In  the  "stationary"  model,  the  plastic  work  done  in  a  single  cycle  of 
amplitude  s  >  Oq  is 

p  2(s-cu)  P  p 

W  =  - p - -V[(E  +aEe)o  +  Es]  ,  (32) 

s  (2Eh+  :e}k  u 
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where  the  notation  is  as  it  was  in  Sec.  2.1.  Since  flow  stress  oq  will  be 

D 

shown  to  change  with  fatigue,  it  is  interesting  to  consider  W  as  a  function 
of  oq,  as  shown  in  Fig.  2.9.  Also  shown  in  Fig.  2.9  is  the  same  function  in 
the  "kinematic"  model,  where 


,,P  _  2o0'S  ’  °0) 
*  =  EP  <■  at6 


(33) 


P  P 

Note  that  W  and  W  always  have  the  same  maximum  value,  although  it  generally 

^  k  p 

occurs  at  different  values  of  the  flow  stress.  When  t  =  0  (the  perfectly 
P  P 

plastic  limit),  W?  and  Wk  are  identical.  Note  that  if  the  flow  stress  we^e  to 
fall  from  a  high  value  near  the  applied  load  to  nearly  zero  during  fatigue, 
the  rate  of  accumulation  of  plastic  work  will  rise  and  then  fal 1 .  This  is 
just  another  manifestation  of  the  opening  and  closing  of  the  hysteresis  loop 
under  the  same  conditions  noted  in  Sec.  2.1. 


teat  3;*4. 
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Fig.  2.9  The  plastic  work  done  in  a  homogeneous  grain  in  one  cycle  during 
equilibrium  hysteresis  in  the  "stationary"  and  "kinematic"  models 
as  a  function  of  flow  stress  on. 
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From  the  interpretation  of  data,  it  can  also  be  conjectured  that 

P  p  D 

plastic  stiffness  E  changes  with  fatigue.  Therefore  W  and  W,  are  plotted  in 

_  O  fN 

Fig.  2.10  as  a  function  of  E  .  There  are  potentially  important  qualitative 

differences  between  the  two  models.  In  the  "kinematic"  model  a  simple  inverse 

relationship  exists.  In  the  “stationary"  model,  the  shape  of  the  function 

depends  on  the  values  of  the  applied  cycle  stress  amplitude  s  and  flow  stress 

P  P 

cq.  If  s  =  2oq,  the  W$  function  is  different  from  only  by  a  numerical 
factor.  However,  if  2og  <  s  <  3cg,  a  point  of  inflection  is  present  in  , 
and  if  s  >  3ag,  a  maximum  exists  at  EP  >  0. 

The  conclusion  to  be  drawn  from  this  brief  comparison  is  that  even 
during  constant  amplitude  loading,  if  the  plastic  stiffness  is  not  constant, 
it  is  important  (for  the  purpose  of  calculating  plastic  work  done)  to  estab¬ 
lish  a  good  model  for  the  yield  criterion. 


(s  -  c„r 


S  >  3 
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Fig.  2.10  The  plastic  work  done  in  a  homogeneous  grain  in  one  cycle  during 
equilibrium  hysteresis  in  the  "stationary"  (solid  lines)  and 
"kinematic"  models  (dashes)  as  a  function  of  plastic  stiffness  or 
work-hardeni ng  parameter. 
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The  addition  of  a  plastic  work  calculation  based  on  Eq.  (31)  to  the 
computer  program  for  the  biaxial  case  is  trivial.  Because  of  the  structure  of 
the  program,  Wp  would  then  be  calculable  for  any  assumed  loading  sequence. 

One  possible  means  of  testing  a  model  for  yielding  is  then  to  compare  computer 
simulations  for  different  loading  sequences  to  see  if  the  predicted  accumula¬ 
tion  of  plastic  work  can  be  correlated  with  experimental  observations  (e.g., 
initiation  events)  that  are  expected  to  be  sensitive  to  plastic  work. 

2 .5  Other  Associated  Problems  Studied 

During  the  life  of  this  program,  several  substantial  tasks  were 
undertaken  that  have  not  been  presented  as  part  of  the  deformation  models  in 
the  preceding  sections.  Nevertheless,  the  information  gleaned  from  these 
various  studies  remains  an  essential  part  of  the  rationale  of  our  current 
models.  Furthermore,  some  of  the  approaches  that  we  discontinued  (usually 
because  of  time  limitations)  remain  viable  and  interesting  ideas  for  future 
work.  We  shall  therefore  discuss  them  here  briefly  before  moving  on  to  report 
the  experimental  work  and  its  analysis. 

2.5.1  The  Linear  Chain  Model 

At  the  beginning  of  the  program,  the  extent  of  interaction  between 
two  neighboring  or  close  deforming  grains  was  unknown.  To  examine  this  ques¬ 
tion,  a  crude  linear  model  of  interacting  grains  was  formulated.  A  linear 
chain  was  chosen  to  represent  the  two-dimensional  surface  of  the  specimen  in 
the  belief  that  the  important  mutual  interactions  under  uniaxial  loading  were 
felt  only  between  grains  lying  in  close  proximity. 

The  mechanism  for  interaction  between  different  grains  was  assumed  to 
be  the  change  in  local  stress  in  one  grain  caused  by  the  elastic  constraint 
field  generated  by  the  deformation  in  the  other.  This  field  was  calculated  by 
assuming  that  each  grain  was  ellipsoidal  and  homogeneously  deformed,  and  by 
invoking  Eshelby's7  solution  for  the  field  outside  an  ellipsoidal  inclusion. 
Coding  this  solution  on  a  computer  constituted  the  major  part  of  the  work  on 
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this  model.  The  change  in  local  stress  in  a  grain  because  of  reaction  to  its 
own  deformation  was,  as  in  the  extant  model  described  in  Secs.  2.1  to  2.4, 
calculated  from  Eshelby’s4  simpler  solution  for  the  field  inside  an  inclusion. 

The  outcome  of  the  chain  model  was  moderately  surprising  but  now 
believed  by  us,  in  the  light  of  much  more  experimental  data,  to  be  clearly  the 
case.  The  model  showed  that,  as  long  as  the  plastic  deformation  in  any  grain 
is  volume  conserving,  which  is  assuredly  the  case  to  a  high  degree  of  preci¬ 
sion,  then  the  change  in  local  stress  in  even  a  neighboring  grain  is  rather 
smal 1 .  In  fact,  compared  to  the  reaction  to  the  grain's  own  deformation,  it 
is  essentially  always  neglible.  Therefore,  it  became  clear  that  most  of  the 
experimental  data  could  and  ought  to  be  interpreted  by  assuming  that  each 
grain  behaves  independently  of  its  neighbors,  as  though  embedded  in  an  en¬ 
tirely  elastic  matrix.  This  is  a  key  premise  of  all  of  the  models  presented 
above.  The  coding  of  Eshelby's  solutions  for  deformation  outside  the  plastic 
grains  later  proved  to  be  essential  in  the  experimental  verification  of  this 
hypothesis  based  upon  a  comparison  of  measured  and  predicted  strains  just 
outside  of  a  deformed  grain. 

It  should  be  stated,  however,  that  there  must  be  circumstances  in 
which  the  model  breaks  down  for  an  isolated  deforming  grain.  If  the  coherence 
of  the  grain  boundary  fails,  or  if  slip  bands  penetrate  into  a  neighboring 
grain,  or  if  plastic  deformation  begins  to  occur  equally  over  groups  of  grains 
or  even  the  whole  surface,  then  a  new  model  is  needed.  In  such  cases,  a  more 
realistic  model  than  the  crude  linear  chain  model  would  also  be  required. 

2.5.2  Corrections  Due  to  the  Free  Surface 

The  reason  that  so  many  interesting  results  are  so  readily  obtained 
from  the  models  described  above  for  deformation  in  a  single  grain,  is  essen¬ 
tially  that  great  simplicity  is  afforded  by  being  able  to  assume  consistently 
that  the  deformation  remains  homogeneous.  However,  this  is  strictly  correct 
only  if  the  grain  is  an  entire  ellipsoid  embedded  in  an  infinite  matrix.  In 
our  case,  we  have  been  modelling  approximately  semi-ellipsoidal  grains  at  a 
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free  surface  of  the  specimen.  It  was  therefore  essential  to  examine  closely 
the  effects  of  the  presence  of  the  free  surface,  and  to  this  end  the  following 
numerical  work  was  carried  out. 


It  was  assumed  that  the  deformation  in  the  grain  consists  solely  of 
pure  extensions  and  pure  contractions  in  the  axis  system  defined  by  the 
grain's  own  axes.  A  free  surface  in  a  plane  of  symmetry  of  the  grain  can  then 
be  created  by  cancelling  the  normal  tractions  found  across  that  plane  in  the 
solution  of  the  constrained  entire  ellipsoid  in  an  infinite  space  (see  Fig. 
2.11).  The  changes  to  the  stress  field  at  any  point  on  the  surface  (or  below 
it,  if  required)  due  to  this  cancellation  can  be  found  by  integrating  Mind- 
lin's^  Green's  function  for  the  half-space.  This  procedure  was  followed 
numerically  by  writing  a  computer  program. 

To  calculate  changes  in  the  in-plane  stresses  o  and  o  (where  the 

*  y 

free  surface  is  the  plane  z  =  0)  ,  one  must  evaluate  the  integral 
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(34) 
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Fig.  2.11  The  cancellation  of  normal  tractions  at  the  surface  z  =  0  to 
form  a  free  surface . 
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where  cz  is  the  normal  stress  to  be  cancelled  at  the  point  (x,y)  and  the  ori¬ 
gin  lies  at  the  point  where  the  correction  is  required.  Quantity  v  is  Pois¬ 
son’s  ratio.  The  singular  nature  of  Eq.  (34),  together  with  the  fact  that  oz 
generally  suffers  large,  rapid  changes  just  outside  the  ellipsoid,  make  the 
numerical  task  quite  difficult.  Nevertheless,  in  the  short  time  available, 
sufficiently  accurate  solutions  were  obtained  to  reach  some  important 
conclusions . 

The  changes  in  the  components  of  the  in-plane  surface  stresses  accom¬ 
panying  the  creation  of  the  free  surface  show  a  singularity  at  the  grain 
boundary.  This  is  the  dominant  feature  of  the  example  shown  in  Fig.  2.12. 
Analytical  work  on  the  surface  displacements  under  another  contract  has  proven 
that  the  accompanying  singularity  in  the  derivative  of  the  normal  displacement 
is  logarithmic,  and  the  implied  weakness  of  the  divergence  is  borne  out  by  the 
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Fig.  2.12  An  example  of  the  change  in  cx  along  the  applied  stress  axis 
due  to  the  creation  of  the  free  surface.  For  a  volume 
conserving  deformation  in  the  uniaxial  model,  normalized  to 
the  value  of  ox  inside  the  entire  ellipsoid  before  the  creation 
of  the  free  surface.  Note  a  certain  amount  of  numerical  noise, 
a  =  203  ^m,  b  =  103  un  and  c  =  50  urn.  There  is  a  logarithmic 
singularity  at  the  boundary  which  is  averaged  in  the  numerical 
procedure.  This  will  be  relaxed  by  localized  plasticity. 
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numerical  results,  which  show  that  it  is  very  limited  in  range.  The  other 
feature  of  Fig.  2.12  is  the  magnitude  of  the  stress  corrections,  which  have 
been  normalized  against  the  internally  constrained  stress  calculated  for  the 
entire  homogeneous  ellipsoid.  They  are  very  small,  even  quite  close  to  the 
grain  boundary.  This  has  always  been  found  to  be  the  case  if  the  plastic 
deformation  in  the  ellipsoid  was  volume  conserving.  If  it  was  not  volume 
conserving,  then  the  corrections  were  found  to  be  quite  substantial,  typically 
ranging  up  to  0.3  on  the  normalized  scale. 

Since  plastic  deformations  are  indeed  volume  conserving,  we  proceeded 
to  ignore  the  free  surface  effects  in  our  modelling.  However,  an  accurate 
analysis  of  strains,  especially  in  the  important  grain  boundary  region,  will 
probably  require  their  inclusion.  We  are  now  in  a  good  position  to  pursue 
this  aspect  in  the  future,  because  of  the  development  of  the  numerical 
sol uti ons . 

2.5.3  Micromechanical  Considerations 

There  are  two  levels  on  which  the  deformation  measurements  can  be 
studied,  differing  in  the  physical  scale  at  which  attention  is  directed.  The 
models  presented  above  belong  to  the  "macroscopic"  scale  in  the  sense  that  we 
have  treated  the  material  as  possessing  idealized,  smoothly  varying  constitu¬ 
tive  properties  without  consideration  of  the  processes  on  a  much  finer, 
"microscopic"  scale  that  lead  to  them.  On  the  microscopic  scale,  in  this 
context,  are  dislocations,  slip  bands,  dislocation  cells,  and  pile-ups.  A 
feasible  alternative  approach  to  interpreting  the  experiments  is  to  start  on 
the  microscopic  scale  in  the  first  place,  building  up  to  the  calculation  of 
measurable  macroscopic  strains. 

This  kind  of  approach  has  been  used  in  the  past  by  many  other 
workers.  One  application  of  particular  interest  is  that  of  Tanaka  and  Mura*1 
in  their  calculation  of  ratcheting  plastic  strains  in  a  deforming  grain.  In 
this  work,  Muskhel  i shvi  1  i  * s  solution^  for  the  singular  integral  equation 
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governing  the  dislocation  density  in  a  one-dimensional  pile-up  is  used.  In 
such  an  equation,  the  dislocations  experience  mutual  interactions,  usually 
taken  from  the  elastic  solution  to  the  field  around  a  dislocation  line  in  an 
infinite  medium.  These  interactions  sum  to  produce  a  back  stress,  which  cor¬ 
responds,  on  the  macroscopic  scale,  to  the  reaction  stress  in  a  constrained 
inclusion  calculated  by  Eshelby's  technique  (as  described,  for  example,  in 
Sec.  2.1).  In  either  case,  the  back  stress  or  reaction  stress  acts  against 
further  plastic  deformation. 

We  have  investigated  the  feasibility  of  solving  a  two-dimensional 
system  of  interacting  dislocations  to  model  deformation  in  a  grain.  A  singu¬ 
lar  integral  equation  in  two-dimensions  was  written  with  piece-wise  continuous 
boundary  conditions  at  the  grain  boundary  to  distinguish  segments  that  were 
intact  (diverging  dislocation  density)  and  segments  where  the  boundary  had 
failed  and  separated  (vanishing  dislocation  density).  Analytical  solutions  to 
such  a  problem  are  not  known,  but  standard  numerical  techniques  have  been  used 
on  analogous  problems  in  fluid  mechanics  with  great  success.  Therefore,  this 
approach  should  be  considered  viable,  and  the  solution  obtained  would  form  a 
very  valuable  link  between  micromechanical  concepts  and  the  phenomenological 
constitutive  relations  used  in  studies  on  the  macroscopic  scale.  We  spent 
some  time  considering  our  options  for  implementation  of  such  a  model  because 
success  would  allow  the  representation  of  deformation  associated  with  grain 
boundary  triple  points,  the  interaction  with  internal  grain  deformation  and 
better  insight  into  the  progression  of  fatigue-induced  changes  in  flow 
stress.  We  concluded  that  while  the  opportunity  was  there,  we  did  not  have 
time  to  pursue  this  line  further  in  this  program.  To  take  it  to  a  useful 
conclusion  would  have  been  a  major  project  in  itself. 
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3.0  MATERIALS  AND  EXPERIMENTAL  PROCEDURES 


3.1  Materi a] s 

The  A1  2219-T851  alloy  has  been  used  extensively  in  prior  studies  of 
crack  initiation  and  short  crack  growth.  Pertinent  properties  of  the  material 
are  its  cyclic  yield  strength  of  360  MPa  and  mean  longitudinal  and  short 
transverse  grain  diameters  of  60  and  20  pm,  respectively.  The  material  has 
the  highest  as-received  hydrogen  content  that  we  have  encountered,  and  is 
therefore  more  microplastic  than  the  average  2219  alloy.  Specimens  were 
fatigued  in  a  dry  nitrogen  atmosphere  to  further  magnify  the  microplastic 
strains  achieved  (see  Ref.  13).  Miniature  tapered  cantilever  beam  specimens 
were  used^  and  were  fatigued  in  fully  reversed  bending  at  +85%  of  the  yield 
strength.  Prior  to  fatigue,  the  surfaces  were  prepared  by  a  careful  machining 
and  polishing  schedule  involving  use  of  progressively  finer  cuts  then  progres¬ 
sively  finer  polishing  grits,  ending  with  an  0.03  pm  MgO  powder.  The  residual 
surface  stresses  left  by  this  procedure  and  measured  by  x-ray  diffraction  are 
small,  +4  MPa.  After  fatigue,  localized  strain  measurements  were  made  vs 
applied  surface  stress  for  values  less  than  or  equal  to  the  maximum  stress 
incurred  during  fatigue,  using  a  bending  jig  in  our  SEM  (see  Ref.  13). 

3 .2  Measurements 

The  strain  measurement  procedures  are  described  el  sewhere , ^ ^  but 
their  essential  features  follow.  The  major  impediment  to  accurate  strain 
measurements  using  a  SEM  is  its  magnification  stability,  which  is  2%  in  our 

O 

instrument.  This  constraint  would  give  a  +2  *  10  strain  measurement  sensi¬ 
tivity  limit  and  is  improved  by  placing  a  floating  reference  or  microscopic 
ruler  on  the  surface  and  gauging  strains  in  the  surface  relative  to  the  ruler. 
The  measurement  involves  comparing  micrographs  of  the  ruler  edges  relative  to 
the  surface  taken  before  and  after  deformation.  We  use  thin  mica  flakes, 
which  are  held  electrostatically  to  the  surface.  Prior  to  fatigue,  these  are 
placed  into  position  within  a  grain  using  a  micro-mani pulator .  Mica  flakes 
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can  be  prepared  which  lie  flat  on  the  substrate  and  are  sufficiently  thin  that 
the  flake  and  substrate  are  co-focal  at  magnifications  of  3  *  10^.  To  avoid 
small  parallax  errors  from  substrate  tilting  under  bending,  all  measurements 
in  a  particular  sequence  are  taken  at  the  same  applied  surface  stress.  Dis¬ 
placement  data  are  obtained  from  micrographs  of  the  before  and  after  deforma¬ 
tion  conditions  using  a  stereoimaging  techni que Before  and  after  defor¬ 
mation  pairs  are  viewed  stereoscopically,  and  small  in-plane  relative  dis¬ 
placements  appear  as  height  differences  between  the  mica  gauge  and  the  sur¬ 
face.  A  photometric  flying  spot  is  used  to  convert  the  three-dimensional 
image  into  linear  displacement  data.  We  chose  this  rather  complex  data  reduc¬ 
tion  technique  because  of  its  sensitivity.  The  human  eye  processes  an  entire 
segment  of  the  flake  edge  image  while  making  the  displacement  determination. 
This  greatly  reduces  statistical  errors.  Furthermore,  because  the  comparisons 
involve  only  the  translation  of  an  edge,  which  is  nearly  identical  in  both 
micrographs,  the  displacement  measurement  sensitivity  is  better  than  the  SEM 
point-to-point  resolution.  By  keeping  surface  contamination  low  through  turbo 
pumping  our  SEM,  differential  displacement  accuracies  of  +30A  have  been  reg¬ 
ularly  achieved.  In  analysis  of  the  displacement  data  to  obtain  strains,  it 
can  be  shown  that  only  the  differential  components  of  the  x  or  y  displacement 
need  be  considered  in  calculating  the  principal  strains  provided  that  local 
rigid  body  rotations  vhich  accompany  the  deformation  are  less  that  1°.  No 
examples  where  this  limit  have  been  exceeded  were  found  in  our  study. 
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4.0  ANALYSIS  OF  EXPERIMENTAL  DATA 

Our  experiments  all  involve  the  fatigue  of  A1  2219-T851,  at  stress 
amplitudes  less  than  or  equal  to  +85%  of  the  alloy  bulk  cyclic  yield  strength. 
Regions  of  substantial  plastic  deformation  of  the  surface  are  found  only  in 
grains  larger  than  120  un.  So  the  plastically  active  grains  have  at  least 
double  the  average  grain  size.  Apparently,  the  early  evolution  of  microplas¬ 
ticity  with  fatigue  is  confined  within  these  grains,  which  act  as  isolated 
regions  of  plasticity  constrained  by  an  external  matrix  of  smaller  and  essen¬ 
tially  elastic  grains.  With  continued  fatigue,  however,  some  of  the  largest 
microplastic  grain  sites  show  deformation  which  extends  beyond  their  bound¬ 
aries  into  the  smaller  neighbors.  This  characteri stic  is  especially  apparent 
at  the  longitudinal  boundary  site,  relative  to  the  stress  axis,  as  indicated 
in  Fig.  2.1.  With  extended  fatigue,  cracks  form  at  such  longitudinal  boundary 
sites,  but  in  2219,  it  is  the  earlier  rupture  of  constitutent  particles  within 
the  grains  which  usually  leads  to  the  crack  initiation  important  to  lifetime. 

As  far  as  possible,  we  have  restricted  our  study  to  grains  whose 
deformation  can  be  analyzed  using  the  models  of  Sec.  2.0.  These  are  grains 
whose  surface  cross-section  is  crudely  elliptical,  with  their  major  axis 
parallel  to  the  loading  axis.  When  looking  at  deformation  at  the  grain  bound¬ 
aries,  we  have  avoided  locations  near  triple  points  and  sought  transverse 
sites  with  boundaries  nearly  parallel  to  the  external  stress  and  longitudinal 
sites  with  boundaries  nearly  perpendicular  to  the  external  stress  axis.  Where 
appropriate  to  the  analysis,  the  assumption  that  the  deformation  within  a 
microplastic  grain  was  uniform  has  been  made  in  obtaining  local  values  of  flow 
stress  from  the  experimental  data. 

This  limitation  appears  to  be  acceptable  as  long  as  the  grain  inter¬ 
ior  is  substantially  softer  than  its  surroundings.  The  dual  domain  deforma¬ 
tion  model  is  applied  to  estimate  flow  stresses  near  the  grain  boundaries. 

Over  small  numbers  of  cycles  the  local  flow  stresses  within  the  plastically 
active  grains  appear  to  be  quasi -static .  Within  the  time  frame  of  a  short 
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loading  segment,  load  sequence  effects  apparently  arise  from  a  coupling 
between  the  mean  stress  in  the  grain  and  the  reaction  stress  to  constraint  of 
its  deformation  aefined  by  the  load  sequence.  Over  longer  fatigue  segments, 
progressive  changes  in  the  local  flow  stress  must  be  considered.  The  fol¬ 
lowing  analysis  deals  with  four'  topics. 

1.  Flow  Stress  at  Grain  Interiors  -  The  flow  stresses  are  examined 
relatively  early  in  fatigue  using  a  uniform  stationary  deforma¬ 
tion  model,  and  applied  to  grains  where  the  deformation  is 
apparently  confined  to  the  grain  interior. 

2.  Flow  Stresses  Near  Grain  Boundaries  -  The  dual  domain  model  is 
used  to  characterize  flow  stress  variations  within  grains  with 
emphasis  on  differences  between  sites  in  the  interior  and  near 
the  boundaries.  [These  results  utilize  NSF  data,  but  the  flow 
stresses  come  from  numerical  calculation  done  by  the  NADC/NAVAIR 
code  described  in  Sec.  2.3.] 

3.  Flow  Stress  Evolution  -  This  subsection  presents  what  is  now 
known  about  the  progressive  changes  in  local  flow  stress  with 
fatigue,  including  the  plastic  disruption  of  grain  boundaries. 

4.  Load  Sequence  Effects  -  The  response  of  local  strain  ranges 
within  microplastic  grains  to  sudden  changes  in  external  load 
are  examined. 

4 .1  Flow  Stresses  at  Grain  Interiors 

First  we  examine  some  properties  of  the  uniform  deformation  of  an 
ellipsoidal  grain  by  considering  the  stationary  deformation  model.  An  impor¬ 
tant  result  is  that  the  experimentally  observed  width  of  the  stress  (external) 
vs  strain  (local)  hysteresis  loop  in  the  grain  is  much  less  than  the  actual 

P 

total  plastic  strain  z  experienced  over  a  loading  cycle.  The  example  (Fig. 
— -  max 
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4.1)  is  for  fully  reversed  loading  with  s  =  305  MPa.  The  loop  width  WQ  is 
measured  at  zero  external  load  (i.e.,  oa  =  0  )  ,  and  eP  is  the  plastic  strain 
range  encountered  from  minimjn  to  maximum  load  in  the  cycle.  Figure  4.2  makes 
the  meaning  of  these  two  quantities  a  little  clearer.  The  stress  range  in  the 
stress  (local)  vs  strain  (local)  loop  in  Fig.  4.2  is  smaller  because  of  the 
reaction  stress  caused  by  constraint  of  the  grain's  microplasticity.  If 
yielding  in  tension  is  at  a  positive  stress,  then  eP  is  just  the  width  of 

JTI3  X 

the  stress  (local)  vs  strain  (local)  loop  at  a  (local)  =  0.  Notice  in  Fig. 

4.1  that  the  measurable  WQ  reaches  a  maximum.  This  happens  when  yielding 
occurs  at  ca  =  0  .  We  are  interested  in  this  behavior  for  its 


4  x  10-3i 


-tpmax*°  =  0.6) 


tpmax<»  =  0.8). 


Wo(0  =  0.81- 


0  100  200  300  400 

FLOW  STRESS.  MPa 

Fig.  4.1  Numerical  values  of  loop  widths  and  plastic  strain  range  obtained 
with  the  uniaxial  stationary  homogeneous  deformation  model  for 
s  =  305  MPa  and  for  several  values  of  a. 
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Fig.  4.2  Stress-strain  response  as  experienced  inside  a  grain  a  (local) 
and  as  apparent  to  an  experimenter  ca  (applied). 

experimental  implications  and  also  because  of  its  pertinence  to  crack,  initia¬ 
tion  models.  If  hysteretic  energy  is  the  driving  force  for  initiation,  then 
W  is  the  important  loop  parameter,  while  a  plastic  damage  model  would  suggest 

description  by  e  ,  . 

J  max 

The  example  in  Fig.  4.1  has  been  constructed  for  an  elastic-perfectly 
plastic  material  (i.e.,  EP  =  0).  If  we  assume  that  yielding  is  controlled  by 
a  von  Mi ses  criterion,  then  p  =  q  =  -1/2  and  Fig.  2.7  tells  us  that  the  common 
range  of  a  values  for  our  A1  2219  alloy  will  be  0.75-0.95,  with  a  =  0.8  being 
the  most  typical  value.  The  reaction  stress  parameter  a  depends  upon  grain 
shape  and  we  arrived  at  these  numbers  by  assuming  that  the  grain  depth  is 
quite  shallow.  Typical  grain  depths  for  the  material  under  stuuy  are  10  to 
20  pm.  For  such  values  of  a  the  maximum  observable  loop  widths  for  a  micro- 
plastic  grain  in  an  elastic  aluminum  matrix  at  ca  =  +  305  MPa  should  be 
WQ  =  1  x  10'3.  Indeed,  this  is  approx imately  the  experimentally  observed 

D 

maximum,  suggesting  that  E  during  microplastic  deformation  may  be  quite 
srnal  1  . 
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Clearly,  however,  we  require  a  more  direct  determination  of  a  and  E 
to  obtain  the  flow  stress  oq  from  our  experimental  data.  By  a  judicious 
choice  of  load  sequence  experiments,  the  measurements  outlined  in  Fig.  4.3  can 
accomplish  this  analysis  giving  a,  EP  and  jq  for  a  "stationary"  deformation 
model  assumption.  The  measurement  begins  by  fatiguing  the  material  at  a 
constant  stress  amplitude  (s)  to  establish  a  quasi-static  value  of  Oq  in  a 
grain.  Two  sets  of  WQ  loop  widths  are  then  measured  at  oa  =  0  ,  and  involve 
characterization  of  WQ  resulting  from  fatigue  at  +  s' ,  where  s'  x  s.  If  WQ  is 
measured  for  the  first  tensile  loop  to  s'  (case  I  in  Fig.  4.3),  then  the 
stationary  model  predicts 


W 


i 


a  + 

0 


aEe(s  -  a  ) 
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Fig.  4.3  Illustration  of  load  sequence  experimental  used  to  obtain  a, 
Ep  and  Cq,  based  upon  a  "stationary"  model  analysis. 
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If  the  loop  opening  for  fatigue  at  +  s'  is  allowed  to  return  to  equilibrium 
after  the  stress  has  been  dropped  from  s  (Case  II  in  Fig.  4.3),  then 


W0(2) 


o0)(l  -  a) 
+  aE6 


(4.2) 


D 

Relative  to  the  intercepts  02 »  a?  and  emax  defined  in  Fig.  4.3,  a,  E  and 
are  obtained  from  Eqs.  (4.1)  and  (4.2)  and  are  simply: 


oq  =  02 » 


(4.3) 


2( 


°1) 


2(oo  "  +  Sax^ 


( S  -  aQ )  ( 1 


max 


a) 


oT 

2 


(4.4) 

(4.5) 


Also,  for  s'  <  o^,  Wq(1)  =  0.  and  for  s‘  <  o2>  WQ(2)  =  0*  Equations  (4.3  - 
4.5)  must  be  rewritten  if  the  quantity,  s^,  defined  in  2.1.3,  is  less  than 
0,  in  which  case  oj  =  0. 


4.1.1  Analysis  o^  Load  Sequence  Data 

Six  to  eight  values  of  loop  widths  must  be  measured  to  obtain  suffi- 

p 

cient  accuracy  to  calculate  E  ,  a,  and  og  for  a  grain,  using  Eqs.  (4.3)- 
(4.5).  We  have  taken  two  such  data  sets  for  fatigue  at  s  of  approximately 
+  300  MPa;  one  set  is  for  a  330  vn  grain  after  500  cycles  (Fig.  4.4);  the 
other  for  a  120  i/n  grain  after  10^  cycles  (Fig.  4.5).  Both  measurements  were 
taken  near  the  center  of  the  grain  parallel  to  the  x  or  stress  axis  and  over 
gauge  lengths  in  excess  of  100  jin.  Please  notice  that  the  strain  scale  is  an 
order  of  magnitude  more  sensitive  for  the  120  i/n  grain.  The  strains  in  this 
grain  were  tiny.  Results  of  our  analysis  of  these  data  using  Eqs.  (4. 3-4. 5) 
are  given  in  Table  1. 
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Table  1 

Deformation  Properties 


Apparent 

Deformation 

Ep 

°0 

Gy  (0.2%) 

Grain  Size 

Depth 

a 

MPa 

MPa 

MPa 

120  un 

<  2  i/n 

0.95  +0.02 

(9.0  +3.5)  x  104 

100  +20 

280  +60 

330  m 

24  i/n 

0.78 

0 

225  +10 

225  +10 

The  apparent  depths  of  the  deformation  have  been  determined  from  Eshelby  s 
equations  by  using  measured  major  and  minor  grain  axes  at  the  surface  and  the 
calculated  value  of  a.  The  values  of  the  parameters  obtained  for  the  330  ym 
grain  are  quite  consistent  with  other  results  for  0g  and  Ep  (described  later) 
obtained  by  less  elegant  means.  For  grains  in  the  300  i/n  size  range,  fatigue 
apparently  causes  local  softening  to  an  approximately  el astic-perfectly  plas¬ 
tic  state.  Flow  stresses  as  small  as  140  MPa  have  been  found,  and  there  is 
less  reliable  evidence  of  even  smaller  values.  The  estimated  24  ^  depth  of 
deformation  is  consistent  with  the  anticipated  depth  of  the  grain,  but  this 
agreement  still  must  be  tested  by  destructive  analysis. 

The  result  for  the  small  120  in  grain  raises  more  questions  than  it 
answers,  from  a  standpoint  of  the  physics  of  localized  deformation.  Actually, 
for  p  =  q  =  -  1/2,  there  is  no  depth  which  will  give  a  =  0.9,  but  depths  under 
2  in  are  consistent  with  this  a  within  the  range  of  experimental  error.  It 
appears  that  deformation  in  the  120  in  grain  may  be  confined  to  a  very  shallow 
layer  near  its  surface.  A  profiling  of  the  true  depth  of  this  grain  is  also 
planned.  But,  other  data  suggests  that  such  a  shallow  depth  of  the  deforma¬ 
tion  is  real. Were  the  values  encountered  here  for  EP  and  og  common  to  all 
120  in  grains,  WQ  values  of  2-4  x  10"4  would  be  anticipated  for  most  grains  of 
this  size.  Instead  the  1  x  10  loop  width  reported  here  is  an  upper  bound  of 
the  experimentally  observed  values. 
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The  last  entry  in  Table  1  deals  with  the  apparent  0.2%  strain  offset 
yield  strength  of  the  grains.  Remember  with  the  stationary  model  og  is  iden¬ 
tified  as  the  stress  below  which  deformation  is  elastic,  (0.2%)  was  cal¬ 
culated  for  each  grain  by  using  EP  and  og  to  construct  a  a  (local)  vs  ex 
(local)  plot  for  the  grain  interior  such  as  in  Fig.  4.2.  The  0.2%  offset  is 
then  determined  in  the  normal  manner. 


4.1.2 


Analysis  of  Old  Loop  Width  Data 


Several  years  ago,  James  and  Morris*3  measured  WQ  at  constant  s  for 
the  same  2219  alloy  studied  here  vs  grain  size,  and  at  several  intervals  in 
fatigue.  It  is  enlightening  to  look  at  these  results  with  the  stationary 
deformation  model.  To  obtain  an  upper  bound  on  the  flow  stress,  we  assume  for 
this  calculation  Ep  =  0  and  an  average  value  of  a  =  0.8.  Because  cQ  vs  Wg  is 
double  valued  in  the  stationary  model  (see  Fig.  4.1),  there  is  an  ambiguity  in 
this  analysis  at  large  values  of  strain.  In  Fig.  4.6,  we  show  results  for 
fatigue  at  +  =  270  MPa  after  104  cycles  and  present  the  maximum  (4.6a)  and 
minimum  (4.6b)  flow  stress  values  found  by  the  model.  Included  in  both  plots 
is  a  datum  (x)  giving  the  bulk  flow  stress,  which  is  assigned  a  location  at 
the  average  alloy  grain  size.  For  the  solid  circles,  we  have  assumed  that  the 
flow  stresses  were  progressively  falling  as  the  plastic  strain  increased.  The 
open  circles  denote  the  ambiguous  stress  values;  in  a)  we  show  the  effect  if 
Og  had  started  to  rise,  in  b)  if  og  continued  to  fall. 

4 .2  Flow  Stresses  Near  the  Grain  Boundaries 


From  limited  data,  it  appears  that  the  flow  stress  within  isolated 
microplastic  grains  may  be  reasonably  uniform.  Elastic  deformation  outside 
the  grain  is  inferred  from  a  comparison  of  measured  strains  vs  distance  from 
grain  boundaries  to  the  results  of  Eshelby  (Sec.  4.2.1).  Just  inside  micro- 
plastic  grains  a  boundary  region  has  been  discovered  which  is  distinguished  by 
a  reversal  of  the  sign  of  WQ  at  oa  =  0  from  that  in  the  grain  interior.  The 
mechanics  of  the  development  of  this  boundary  region  are  explained  by  an 
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Fig.  4.6  Values  of  local  flow  stress  for  A1  2219-T851  fatigued  104  cycles  at 
+  270  MPa  vs  grain  size.  Calculated  for  an  elastic  perfectly 
plastic  grain  using  the  "stationary"  deformation  model,  a)  An 
upper  bound  to  og,  b)  the  most  probable  cQ  value.  The  curve  in  b) 
is  a  (1/grain  size)  function. 
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analysis  provided  by  the  dual  domain  deformation  model.  The  reversal  in  sign 
of  W0  is  a  consequence  of  a  flow  stress  in  the  domain  close  to  the  boundary 
being  slightly  higher  than  the  grain  interior  (Sec.  4.2.2).  Measurements  show 
that  outside  the  microplastic  grain,  the  local  flow  stress  rises  abruptly  to 
that  consistent  with  the  bulk  value. 


4.2.1  Elastic  Matrix  Response 

Except  for  highly  deformed  sites  discussed  in  Sec.  4.3,  the  deforma¬ 
tion  outside  the  isolated  microplastic  surface  grains  is  indistinguishable 
from  elastic.  Reaction  strains  resulting  from  the  confine,  plasticity  have 
been  compared  to  those  calculated  using  Eshelby's  equations  for  several 
grains,  of  which  results  in  Fig.  4.7  are  an  example.  Figure  4.7a  is  of  ex 
along  the  transverse  grain  axis  and  4.7b  gives  along  the  longitudinal 
grain  axis  (see  Fig.  2.1).  The  solid  curves  in  the  elastic  region  just  out¬ 
side  the  grains  have  been  calculated  from  the  strains  just  inside  the  grain 


DISTANCE  FROM  y  BOUNDARY,  Mm 


DISTANCE  FROM  x  BOUNDARY,  Mm 

Fig.  4.7  Example  values  of  strain  just  outside  microplastic  grains 
indicating  an  elastic  matrix  and  deformation  depths  around 
10  uni  -  the  typical  grain  depth.  NSF1  data. 
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for  several  assumed  depths  of  the  plasticity  and  suggest  that  the  microplastic 
deformation  depth  is  comparable  to  the  grain  depth.  We  estimate  that  the  sen¬ 
sitivity  of  our  strain  measurements  precludes  the  distinction  between  values 
of  oq  outside  the  grain  which  are  larger  than  =  300  MPa. 

A  completely  unanticipated  observation  is  that  strains  in  the  micro¬ 
plastic  grain  at  the  transverse  boundary  (Fig.  4.7a)  are  compressive  after 
tensile  loading.  The  results  shown  are  after  a  half  tensile  loading  cycle. 
Tensile  strains  would  be  anticipated  from  a  uniform  stationary  deformation 
model.  About  30  additional  such  measurements  made  under  NSF  funding  confirm 
this  result.  In  the  grain  interior,  WQ  is  positive,  but  a  boundary  layer  some 
20  pm  in  width  is  found  in  which  WQ  is  negative.  An  explanation  of  this  beha¬ 
vior  is  found  in  analysis  of  stress-strain  response  vs  applied  stress  using 
the  dual  domain  model . 

4.2.2  Flow  Stress  vs  Location  Within  a  Grain 

By  means  that  are  explained  momentarily,  we  have  obained  flow  stres¬ 
ses  (og)  shown  in  Fig.  4.8  along  the  y  axis  of  a  270  pm  grain  after  fatigue 
cycles.  In  setting  the  dashed  line,  we  allowed  ourselves  to  believe  that  the 
data  are  consistent  with  a  og  equal  to  the  bulk  flow  stress  just  outside  the 
microplastic  grain.  Also,  since  the  value  of  E**  =  1  x  10+^  MPa  is  found 
inside  the  grain,  ag  -  (0.2%).  The  boundary  layer  mentioned  in  Sec.  4.2.1 

is  only  apparent  in  the  strain  loop  width  (WQ)  (see  Fig.  4.9).  We  have  only  a 
few  points  and  have  used  our  experience  from  strain  measurements  in  other 
areas  to  place  a  dotted  line  through  the  values  of  WQ  measured  at  oa  =  0  after 
a  tensile  loading  cycle. 

An  explanation  for  the  reversal  in  sign  in  W0  at  the  boundary  is 
found  in  the  analysis  of  loop  widths  measured  at  sites  in  the  grain  vs  oa . 

The  experiment  involved  the  measurement  of  strains  over  a  partial  tensile 
cycle,  beginning  and  ending  at  the  same  value  of  ca  (Fig.  4.10).  The  reversal 
in  sign  of  WQ  at  the  boundary  occurs  automatically  in  the  dual  domain  model  if 
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Fig.  4.8  Values  of  og  extracted  from  local  loop  width  data  for  a 
270  urn  grain  using  the  "dual  domain"  model. 
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Fig.  4.9  Illustration  of  the  location  of  the  strain  anomaly  at  a  grain 
boundary  and  flow  stress  values  associated  with  the  anomaly. 
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EXTERNAL  STRESS.  MPa 

Fig.  4.10  Experimental  values  of  loop  width  vs  the  external  stress 

at  the  beginning  and  completion  of  the  loading  cycle.  Data 
are  for  three  different  locations  and  curves  are  from  the 
"dual  domain"  model,  with  E  =  1  x  103  MPa  and  the  flow 
stress  indicated.  NSF1  data. 

oq  at  the  boundary  is  greater  than  in  the  interior  of  the  grain,  provided  that 
EP  is  within  a  very  narrow  range  of  values  near  1  x  io^  MPa,  and  also  that  a 
at  the  boundary  is  small.  Under  this  condition,  the  stress  (external)  - 
strain  (local)  hysteresis  loop  can  have  the  complex  shape  shown  in  Fig.  4.11. 
This  is  the  loop  shape  which  produces  the  result  in  Fig.  4.10b.  At  oa  =  0, 
the  motion  about  the  loop  is  anticlockwise,  but  at  large  joa|  the  motion  is 

p 

conventional.  An  extensive  search  of  the  possible  combinations  of  E  ,  a  and 
oq,  which  give  quantitative  agreement  with  the  data  of  4.10  has  hr .n  used  to 
establish  error  bars  on  the  values  of  Oq  given  in  Fig.  4.8.  Anomalously  small 
values  of  a  are  needed  to  describe  the  experimental  results,  however,  and  sug¬ 
gest  that  the  boundary  constraint  is  still  not  fully  understood.  We  consider 
this  aspect  later  (Sec.  5.0). 
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Fig.  4.11  External  stress-local  strain  behavior  at  a  boundary  harder 
than  a  grain  interior  for  appropriate  a  and  E  ,  predicted 
by  the  dual  "domain  model"  has  a  characteristic  structure 
useful  for  analysis  of  local  flow  stress  values. 

4 .3  Flow  Stress  Evolution 

We  do  not  yet  have  a  model  for  the  variation  of  local  flow  stress 
with  fatigue,  eg  must  be  related  both  to  grain  size  and  to  load  spectra. 
However,  an  analysis  (Fig.  4.12)  of  localized  strain  data  obtained  at  constant 
amplitude  does  provide  some  insight  into  trends  in  og  with  grain  size  and 
fatigue  cycles  in  the  A1  2219-T851  alloy.  The  values  given  in  Fig.  4.12  have 
been  extracted  from  strain  measurements  previously  obtained  by  James  and 
Morris,*3  by  using  a  homogeneous  stationary  deformation  model.  The  analysis 
involves  several  assumptions  which  remain  to  be  tested.  In  particular,  for 
large  grains  the  stress-strain  loop  widths  are  observed  to  increase  to  a 
maximum  and  then  decrease  with  progressive  fatigue.  For  the  350  urn  grain  case 
in  Fig.  4.12,  we  have  taken  the  flow  stress  consistent  with  observed  loop 
width  values,  by  assuming  that  og  continues  to  fall  after  the  maximum  WQ  is 
reached.  However,  there  are  three  other  possibilities,  namely 
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Fig.  4.12  Calculated  values  of  local  flow  stress  with  fatigue  obtained  by 
averaging  values  for  a  number  of  grains  of  sizes  near  those 
indicated.  Values  less  than  about  140  MPa  are  less  credible 
because  they  are  very  model  sensitive. 

a)  oq  itself  reaches  a  minimum  and  then  increases; 

b)  Ep  increases; 

c)  The  grain  boundaries  soften,  decreasing  the  applied  stress  range 
in  the  grain  center. 

We  have  looked  at  several  special  cases  of  very  large  grains  (330  pm 
and  larger)  and  it  is  clear  for  these  that  plastic  deformation  of  the  bound¬ 
aries  and  the  neighboring  grains  is  partially  responsible  for  the  observed 
loop  closings.  The  breakdown  of  the  boundaries  is  accompanied  by  violent 
variations  in  localized  strains  which  are  associated  with  slip  banding  in  the 
2219  (Fig.  4.13)  and,  for  very  large  grains  (~  500  pm),  WQ  in  the  grain 
interior  measured  over  a  tensile  loading  cycle,  actually  becomes  compressive 
(Fig.  4.14).  This  is  behavior  reminiscent  of  that  in  the  hard  boundary  zones 
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Fig.  4.13  Residual  component  of  e  induced  by  tensile  loading  after  fatigue 
has  begun  to  disrupt  a  grain  boundary.  NSF1  data. 
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Fig.  4.14  Progress  fatigue  indurad  change  in  WQ  for  two  large  grains  for 
s  =  +  305  MPa.  The  compressive  response  after  500  cycles  shows 
that  the  deformation  is  highly  nonuniform,  and  that  the  grain 
interior  is  probably  harder  than  the  surrounding  matrix. 
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described  previously  and  suggests  that  a  condition  has  been  reached  in  which 
the  grain  interior  has  hardened  relative  to  the  surrounding  medium. 

4 .4  Response  to  Spectrum  Loads 

The  bulk  of  our  information  on  the  effect  of  spectrum  loads  has 
already  been  shown  in  Sec.  4.1.  Initially,  the  only  thought  behind  the  data 
collection  was  to  provide  a  basis  for  testing  of  the  local  deformation  models 
vs  load  sequence.  Later,  we  realized  that  with  an  appropriate  analysis,  the 
response  to  step  changes  in  load  could  also  be  used  to  acquire  estimates  of 
the  important  parameters  representing  the  deformation  properties  of  the  sur¬ 
face.  At  the  same  time,  it  is  clear  from  the  concurrence  of  the  results  with 
theory  that  major  features  of  local  stress-strain  behavior  are  influenced  in  a 
predictable  way  by  local  reaction  stresses  created  in  response  to  the  plastic¬ 
ity.  A  plastic  grain  acts  in  a  fashion  similar  to  a  microscopic  notch. 
Numerical  calculations  have  been  made  of  the  deformation  of  plastic  grains 
under  positive  and  negative  mean  stresses  and  a  ratcheting  of  the  local  defor¬ 
mation  to  achieve  a  local  state  of  fully  reversed  loading  is  found,  just  as 
would  occur  at  a  notch.  This  response  for  surface  deformation  is  consistent 
with  x-ray  diffraction  observations  of  surface  residual  stresses  in  aluminums 
by  James  and  Morris.^  A  next  step  in  testing  understanding  of  this  behavior 
would  be  to  collect  strain  response  data  within  grains  under  nonzero  mean 
loading.  However,  this  research  needs  to  be  carefully  integrated  with  an 
examination  of  the  nature  of  the  stress-strain  and  yielding  behavior  of  the 
fatigue  softened  areas.  As  noted  in  Sec.  2.0,  we  have  to  this  point  relied 
heavily  on  a  "stationary"  yielding  assumption  and  the  analytic  and  experimen¬ 
tal  tools  are  now  in  hand  to  test  this  assumption  and  to  compare  its  predic¬ 
tions  to  those  from  a  "kinematic"  yielding  model.  Remember  that  one  differ¬ 
ence  between  the  two  models  is  that  there  is  no  intrinsic  Baushinger  effect  in 
the  stationary  model.  The  Baushinger  response  seen  when  the  stationary  model 
is  invoked  stems  from  the  reaction  stresses  developed  when  the  plastic  defor¬ 
mation  of  a  grain  is  confined  by  an  elastic  matrix.  On  the  other  hand,  the 
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kinematic  model  material  has  its  own  internal  memory  of  ^ne  stress  history  to 
which  it  has  been  exposed.  Interestingly,  however,  the  strain  in  a  grain 
subject  to  a  kinematic  yielding  assumption  reaches  a  new  equilibrium  in  just 
one  complete  cycle  if  a  step  change  in  the  cyclic  stress  amplituoe  is  made. 
Experimental  observations  (Fig.  4.15)  indicate  that  several  cycles  are  act¬ 
ually  required  for  such  an  equilibrium  to  be  achieved.  However,  before  con¬ 
clusions  are  drawn  from  this  observation,  the  computer  simulations  of  response 
under  the  locally  biaxial  stress  state  caused  by  plastic  deformation  should  be 
examined  carefully  in  the  light  of  additional  experimental  observations. 


Fig.  4.15  Transient  response  in  loop  widths  after  a  step  change  in  cyclic 
stress  amplitude.  In  this  example,  the  initial  stress  range  was 
+  305  Mpa  and  was  then  dropped  to  +  156  MPa,  in  one  experiment, 
and  from  +  305  MPa  to  +  235  MPa  in  a  second  experiment. 
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5.0  DISCUSSION 

We  believe  that  we  have  now  identified  the  major  building  blocks 
which  must  be  included  in  a  model  of  crack  initiation  in  fatigue  under 
spectrum  loads,  at  least  for  structural  aluminum  alloys.  One  key  is  that  the 
softening  of  the  surface  which  occurs  with  fatigue  is  extreme  and  localized  to 
large  grains.  The  values  we  obtain  for  local  flow  stresses  are  to  a  degree 
dependent  upon  the  deformation  model  used  to  analyze  the  data,  but  there  is  no 
doubt  that  the  stresses  become  very  small.  For  A1  2219-T851  values  of  oQ  be¬ 
low  140  MPa  are  credible,  whereas  the  bulk  cyclic  yield  strength  is  360  MPa. 
Furthermore,  the  material  in  a  softened  grain  responds  in  nearly  an  elastic- 
perfectly  plastic  fashion.  Seeing  this,  it  is  little  wonder  that  no  endurance 
limit  has  been  found  for  aluminum  alloys.  Only  under  very  severe  conditions 
such  as  for  large  cyclic  stresses  (perhaps  85T  of  the  bulk  elastic  limit),  and 
large  grains  (perhaps  350  }im) ,  is  there  clear  evidence  that  plasticity  extends 
beyond  a  microplastic  grain  into  neighboring  grains.  In  addition,  it  is  found 
that  the  local  flow  stresses  change  slowly  with  fatigue. 

Over  the  years,  we  have  had  a  chance  to  look  in  detail  at  the  crack 
initiation  behavior  in  aluminum  alloys  and  especially  in  the  A1  2219-T351 
material  used  in  this  study.  For  constant  amplitude  fatigue,  at  lifes  in 
excess  of  10^  cycles,  crack  initiation  in  the  2219  is  entirely  at  constituent 
particles  which  lie  within  large  grains  on  a  specimen  surface.  Only  at  short 
lifetimes  is  grain  boundary  cracking  apparent,  and  normally  it  is  still  the 
constituent  particle  cracking  which  produces  the  fatal  crack.  To  us  this 
suggests  that  the  latter  stage  of  deformation  of  the  surface  in  which  the 
boundaries  of  and  neighbors  to  a  large  grain  begin  to  deform  does  not  neces¬ 
sarily  have  a  consequence  to  lifetime.  Its  pertinence  depends  upon  the  avail¬ 
ability  of  competing  sites  for  initiation  in  the  grain  interiors  which  soften 
first.  Only  in  compression  dominated  spectra  have  we  found  boundary  cracking 
in  2219  to  be  of  special  importance. 


NADC-87046-60 


With  just  these  insights,  we  outline  the  structure  of  a  model  of 
crack  initiation  for  spectrum  loading  and  delineate  the  questions  which  must 
still  be  answered  to  complete  and  test  the  model  (5.1).  Following  that,  we 
examine  what  we  have  learned  that  bears  upon  answers  to  these  questions  for  A1 
2219-T851  (5.2). 

5.1  Future  Initiation  Modelling 

Initiation  modelling  requires  a  physically  correct  representation  of 
localized  microplastic  deformation.  Some  approximations  are  unavoidable  and 
should  be  used  if  they  are  supported  by  experimental  evidence.  We  suggest  the 
following  model  of  crack  initiation  inside  an  isolated  microplastic  grain  is 
appropri ate . 


Grain  Properties  and  Deformation  Model 

Grain  embedded  in  an  elastic  matrix 
Ellipsoidal  grain 

Uniform  deformation  -  response  controlled  by  internal 

reaction  stress 


Material  Properties  with  the  Microplastic  Grain 


Approximately  elastic-perfectly  plastic 
Cycle  dependent  flow  stress 


Need 


[Yield  Model] 

[oq  cycle  and  grain  size  dependence] 
[More  deformation  depth  information] 
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Initiation  Criterion  -  Possibilities 

Need  Peak  plastic  strain  range 

Criterion  Integrated  plastic  strain 

Integrated  hysteretic  energy 

i 

If  the  fatigue  induced  softening  of  a  grain  continues  progressively, 
then  the  peak  and  integrated  plastic  strains  will  follow  suit.  But,  the 
hysteretic  energy  per  cycle  will  decrease  with  sufficient  softening  (as  shown 
in  Sec.  2.4).  Our  particular  interest  in  this  latter  observation  is  that  in 
2219  crack  initiation  in  large  grains  tends  to  occur  in  a  sudden  burst  of 
activity,  if  a  specimen  is  fatigued  at  constant  amplitude.  If  cracking  ceases 

D 

because  E  or  og  begins  to  increase  at  some  point  during  fatigue,  then  the 
same  change  in  initiation  behavior  will  be  seen  at  all  amplitudes.  But,  if 
softening  and  hysteretic  energy  are  involved,  this  cessation  of  initiation 
will  pertain  only  to  constant  amplitude.  The  loop  width  equations  for  the 
"stationary"  deformation  model  show  that,  instead,  a  low-high  amplitude  block 
loading  sequence  will  be  much  less  damaging  than  a  high-low  sequence.  This  is 
a  trend  which  is  supported  by  experiment  for  2219. ^  So  an  important  step  in 
the  model  development  is  to  determine  the  driving  force  for  initiation.  This 
can  easily  be  done  once  progressive  changes  in  og,  Ep  and  W0  have  been  deter¬ 
mined  within  selected  grains  during  fatigue,  using  the  technique  described  in 
Sec .  4.1 . 

The  theory  then  requires  a  model  to  calculate  from  cycle-by-cycle 
local  strain  values  the  fatigue  induced  change  in  og,  and  a  correct  represen¬ 
tation  of  the  local  strains  encountered  in  any  arbitrary  loading  sequence. 

The  latter  is  essentially  provided  by  the  procedures  described  in  Sec.  2.3. 
What  is  still  missing  is  knowledge  of  the  yield  criterion.  For  a  material 
which  is  nearly  an  el astic-penfectly  plastic  material,  such  as  we  find  for 
softened  2219  grains,  the  "stationary"  and  "kinematic"  deformation  models 
should  demark  the  two  extremes  of  the  expected  deformation  behavior.  In  fact, 
if  E  =  0,  the  two  yield  criteria  are  identical.  At  this  juncture,  one  needs 
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to  look  more  carefully  at  the  biaxial  strain  predictions  of  these  models  to 
identify  critical  experiments  which  would  be  useful  for  model  validation. 

5 .2  Local  Mechanical  Properties 

The  solution  we  seek  for  prediction  of  fatigue  initiation  lifetime 
has  been  reduced  to  that  of  a  problem  in  micromechanics.  Models  and  their 
numerical  implementations  to  calculate  local  stress-strain  behavior  in  a 
microplastic  grain  have  been  developed,  and  for  fully  reversed  loading,  a 
"stationary"  local  yielding  criterion  provides  an  acceptable  representation  of 
experimental  observation.  Fulfillment  of  our  initial  goal  requires  a  better 
understanding  of  the  way  in  which  the  material  mechanical  properties  (such  as 
flow  stress)  and  the  local  yielding  criterion  evolve  with  fatigue.  A  calcula¬ 
tion  of  the  local  stress-strain  response  in  an  ellipsoidal  grain  for  arbitrary 
loads  awaits  only  this  final  insight.  Although  we  know  that  the  deformation 
within  a  grain  can  be  complex,  especially  near  boundaries,  it  appears  that  an 
assumption  of  uniform  deformation  within  the  grain  interior  will  suffice  for 
our  modeling  objectives  so  long  as  the  important  crack  initiation  is  actually 
at  such  interior  sites. 

As  described  by  Margolin^7  and  Asaro1^  several  years  ago,  any  local 
inhomogeneity  in  deformation  automatically  provides  the  foundation  for  Baush- 
inger  and  load  history  effects.  These  arise  because  of  reaction  stresses 
generated  by  constraint  of  the  nonuni  form  deformation.  For  microplastic 
grains  in  2219  we  find  that  these  constraining  stresses  are  huge  because  the 
local  values  of  flow  stress  become  very  small  during  fatigue.  Herein  is  the 
major  problem  in  analysis  of  surface  microplastic'ty.  The  local  and  applied 
stresses  differ  substantially  and  only  the  strains  and  applied  stress  are 
experimentally  available.  Appropriate  theories,  experimentally  validated,  are 
needed  to  discern  the  local  mechanical  properties  and  yielding  behavior  almost 
hidden  from  view  by  the  reaction  stresses. 

This  is  what  we  now  know  about  the  mechanical  properties  of  micro- 
plastic  grains  in  2219.  The  most  reliable  values  for  properties  in  a  small 
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(120  pm)  grain  are  reported  in  Table  1.  Although  the  data  were  obtained  after 
fatigue,  we  think  it  is  likely  that  the  numbers  represent  the  ambient  prefa¬ 
tigue  values  of  the  surface.  It  was  necessary  to  measure  surface  strains  in 
the  10“^  range  to  obtain  the  parameter  values  of  Oq  =  100  MPa  and  Ep  = 

1 .3  *  104  MPa  for  this  example.  Bulk  materials  are  known  to  have  a  two-step 
yielding  behavior  and  the  100  MPa  value  would  correspond  to  the  lower  yield 
point  for  the  grain.  It  happens  that  the  apparent  depth  of  the  layer  with 
this  flow  stress  is  only  2  pin;  this  is  typical  of  the  dislocation  cell  size  in 
A1  -  suggesting  a  possible  connection.  This  would  mean  that  EP  for  the  cell 
structure  would  be  in  excess  of  104  MPa.  Values  of  og  for  100  pm  grains  given 
in  Fig.  4.6  are  quite  consistent  with  this  result.  Remember,  the  estimates 
shown  in  Fig.  4.6  are  actually  upper  bound  estimates  of  aQ  in  that  we  have 
assumed  that  EP  =  0,  independent  of  grain  size  in  calculating  the  Og's 
given.  Supporting  evidence  that  the  mechanical  properties  found  for  the 
120  pm  grain  may  be  universal  in  2219  is  apparent  in  Fig.  4.4.  In  these  loop 
width  results  for  a  330  pm  grain,  a  small  positive  plastic  strain  is  found 
below  the  s'  =  a2  axis  intersection  point,  contrary  to  the  "stationary"  defor¬ 
mation  model  Wq(2)  =  0  prediction.  It  happens  that  the  very  small  residual 
strain  present  is  just  those  that  would  be  anticipated  if  there  were  an  ear¬ 
lier  yield  point  in  the  grain  with  a  plastic  modulus  of  about  Ep  ~  104  MPa  for 
a  stress  below  . 

The  fatigue-induced  softening  of  the  surface  of  2219  appears  to 
involve  the  development  of  long-range  slip  extending  grain  wide.  The  bands 
are  not  seen  by  optical  microscopy,  but  are  apparent  in  strain  field  measure¬ 
ments  made  over  the  entire  surface  of  a  grain  [NSF*  result].  We  think  it  is 
this  longer  range  dislocation  structure  which  exhibits  an  elastic-perfectly 
plastic  (i.e.,  EP  =  0)  deformation  behavior.  A  recent  paper  by  Laird,  to  be 
published  shortly,  shows  that  the  fatigue  of  dispersion-hardened  single  copper 
crystals  produce  banding  and  a  drop  in  flow  stress  to  a  value  approaching  that 
of  the  unalloyed  copper.  This  happens  for  a  plastic  strain  in  the  range  of 
1-5  x  l o- 3 ,  numbers  typical  of  those  seen  for  A1  in  this  study.  Only  at 
higher  plastic  strains  is  hardening  and  the  development  of  a  tangled  disloca- 


NADC-87046-60 


tion  structure  apparent.  If  the  same  phenomenon  were  active  in  the  single- 
grain  crystals  of  2219  at  the  surface,  it  would  produce  essentially  the  same 
characteristic  changes  in  local  flow  stress  with  fatigue  that  we  have  reported 
here.  Only  the  dynamics  of  the  softening  must  be  more  complex,  because  in  a 
surface  grain  the  local  plastic  strain  range  changes  cycle-by-cycle.  Never¬ 
theless,  the  net  effect  on  the  local  ag  should  be  the  same. 

The  mechanical  properties  near  grain  boundaries  in  2219  are  more 
nebulus.  Clearly,  early  in  fatigue  the  boundaries  are  harder  than  the  grain 
centers.  But,  despite  the  obvious  success  of  the  dual  domain  model  in  simu¬ 
lating  boundary  stress-strain  response,  there  is  an  inconsistency  in  the 
model.  The  problem  concerns  a(2),  the  relative  reaction  stress  parameter  at 
the  boundary.  To  explain  our  observations,  approximate  values  for  a(2)  of  0.7 
are  required.  While  boundary  zone  structures  can  readily  exist  which  produce 
such  an  a(2),  there  is  no  obvious  zone  shape  which  would  create  such  a  layer 
along  the  entire  length  of  the  boundary.  Yet,  that  is  exactly  what  is  ob¬ 
served.1  It  appears  that  application  of  the  free  surface  correction  to  our 
model  does  lower  the  boundary  a  values  somewhat,  but  probably  not  enough  to 
explain  what  we  see.  An  additional  consideration  is  that  the  triple  points  in 
the  boundary  are  also  soft  spots.1  The  fact  is  that  deformation  near  the 
boundaries  is  complex  and  is  still  poorly  understood. 

A  second  key  to  understanding  the  behavior  of  deformation  localized 
in  a  single  grain  is  that  the  measurable  hysteresis  loop  widths  in  the  grain 
are  very  small  because  of  the  constraint  provided  by  the  surrounding  elastic 
matrix.  Despite  this,  the  local  plastic  strains  can  be  large.  They  must 
approach  1  x  io"^  in  some  grains  we  have  examined.  This  means  that  the  elas¬ 
tic  reaction  stresses  in  a  grain  are  also  very  large.  Herein,  apparently,  is 
the  source  of  many  load  amplitude  and  sequence  effects.  The  mean  stress  and 
stress  range  experienced  in  the  elastic-plastic  material  within  a  grain  is 
derived  from  both  the  external  stress  and  from  a  reaction  stress  which  depends 
upon  the  immediate  state  of  plastic  strain  in  the  grain. 
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Much  of  what  was  previously  known  about  the  properties  of  grain  inte¬ 
riors  vs  boundaries  in  metals  comes  from  a  few  monotonic  deformation  studies. 
Except  at  triple  points,  the  boundaries  are  typically  found  to  be  harder  than 
the  grain  centers.  But,  these  are  bulk  deformation  observations  and  the 
applied  plastic  strains  are  larger.  Even  the  differences  in  flow  stress  found 
between  boundaries  and  the  grain  interior  have  been  relatively  small.  So  our 
findings  regarding  the  fatigue-induced  softening  of  A1  2219-T851  are  a  major 
revelation  of  the  physics  and  mechanics  of  surface  deformation.  The  results 
are  essential  to  an  understanding  of  crack  initiation.  They  also  have  immedi¬ 
ate  application  to  the  understanding  of  fatigue-induced  changes  in  residual 
surface  stresses.  Additionally,  they  provide  a  probable  explanation  for  the 
persistence  of  propagation  of  short  cracks  in  a  surface  even  at  cyclic  stress 
intensities  well  below  the  long  crack  stress  intensity  range  threshold.  Now 
we  realize  that  the  surface  is  not  even  approximately  elastic.  The  small  sur¬ 
face  cracks  initiate  in  large  grains  and  are  embedded  in  a  highly  microplastic 
surface  layer.  These  insights  have  come  in  quite  a  rush  as  this  program  has 
drawn  to  completion  and  papers  detailing  our  results  are  still  in  preparation. 
The  first,  entitled,  "Microplastic  Surface  Deformation  of  A1  2219-T851,"  will 
be  submitted  to  Acta.  Met. 
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6.0  SUMMARY  AND  CONCLUSIONS 

Fatigue  cracks  in  structural  aluminum  alloys  commonly  initiate  at  tne 
surface.  Lifetimes  depend  on  the  alloy  microstructure,  especially  on  grain 
size,  and  show  a  distribution  apparently  controlled  by  stochastic  cracking 
processes  which  are  associated  with  individual  grains  at  the  surface.  Conven¬ 
tional  cumulative  damage  methods  are  applicable  to  aluminums,  but  there  are 
cases  where  they  do  not  work  well.  Also,  they  are  especially  inept  in  pre¬ 
dicting  stochastic  behavior  for  spectrum  loading,  and  are  useless  in  a  priori 
predictions  of  the  consequences  of  variability  in  the  material  properties. 

With  the  goal  of  defining  a  methodology  for  making  initiation  lifetime  predic¬ 
tions  on  the  basis  of  physical  insight  into  the  micromechanics  of  the  local 
deformation  processes,  this  research  has  focused  on  the  properties  of  the 
localized  regions  in  the  surface  at  which  cracks  initiate.  This  has  involved 
rather  sophisticated  measurements  of  localized  strains  in  the  surface  in  indi¬ 
vidual  grains,  and  frequently,  over  very  short  gauge  lengths.  But,  the  main 
thrust  of  the  program  has  been  the  development  of  the  theoretical  tools  to 
analyze  our  experimental  data.  The  major  challenge  is  that  the  localized 
plastic  strains  are  large,  and  reaction  stresses  created  by  constraint  of  the 
plasticity  are  also  large.  As  a  consequence,  the  local  stresses  differ  sub¬ 
stantially  from  the  applied,  making  it  impossible  to  directly  measure  local 
yield  strengths  and  plastic  moduli  within  grains  under  study.  The  theoretical 
models  of  deformation  are  needed  to  extract  from  measurements  in  a  grain  of 
total-strain  vs  external  load,  the  true  plastic  strains  and  the  material 
mechanical  properties  which  vary  from  grain  to  grain.  The  models  also  need  to 
establish  from  the  data  the  mechanisms  by  which  these  local  mechanical  proper¬ 
ties  change  with  fatigue  and  also  the  criterion  which  controls  local  yielding 
within  a  grain.  These  analysis  goals  can  only  be  accomplished  by  a  model 
which,  although  it  may  be  an  approximation,  embodies  the  essential  physics  and 
mechanics  of  the  deformation  process.  The  advantage  of  this  approach  is  that 
once  the  mechanical  property  and  yielding  behavior  is  understood,  the  deforma¬ 
tion  model  then  becomes  the  crack  initiation  model.  The  theory  is  used  to 
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calculate  the  local  stress-strain  response  within  a  grain  subject  to  micro- 
plasticity  and  under  conditions  of  arbitrary  load  spectra.  All  that  is 
required  to  complete  the  initiation  lifetime  prediction  is  knowledge  of  the 
crack  initiation  criterion.  Criteria  such  as  critical  values  of  accumulated 
total  plastic  strain,  or  cumulative  local  hysteresis  energy  per  cycle  are 
postulated  and  validated  by  experiment.  Modelling  of  the  stochastic  behavior 
can  then  be  approached  by  numerical  techniques  which  sample  the  lifetime 
response  conditioned  on  a  distribution  of  grain  sizes  and  critical  initiation 
parameten  values.  By  this  point,  however,  the  treatment  of  the  difficult 
spectrum  loading  aspects  of  the  problem  is  already  directly  embedded  in  the 
model  of  localized  deformation  and  requires  no  further  attention. 

We  have  made  substantial  progress  in  the  development  of  such  a 
formalism.  The  first  generation  of  deformation  models  are  complete  and  they 
have  lead  to  startling  discoveries  about  the  mechanical  properties  of  the 
surface  of  the  structural  A1  2219-T851  alloy  which  has  been  the  model  material 
used  in  our  study.  We  have  also  gained  substantial  insight  into  the  condi¬ 
tions  under  which  the  approximations  in  the  deformation  models  will  be  valid, 
and  information  on  the  way  the  localized  mechanical  properties  within  grains 
change  with  fatigue.  Needed  to  complete  the  initiation  methodology  is  a 
physical  model  of  these  property  changes  in  relation  to  local  grain  size  and 
to  the  local  stress-strain  environment.  Also  more  must  be  learned  about  the 
local  yield  criterion,  and  experiments  are  needed  to  establish  which  is  the 
most  appropnate  local  crack  initiation  criterion.  Thereafter,  remaining  work 
on  the  predictive  methodology  will  simply  be  numerical  modelling  to  collate 
these  insights  and  validation  studies  of  the  lifetime  predictions. 

In  summary,  the  following  has  been  accomplished. 
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Theoretical  Modeling 

A  numerical  implementation  of  Eshelby's  equations  to  calculate  the  elastic 
strain  and  stress  fields  surrounding  a  plastically  deformed  grain  has  been 
completed.  Results  for  selected  grain  shape  showed  that  the  fields  exter¬ 
nal  to  the  grains  might  be  sufficiently  large  to  damage  regions  near  the 
boundaries,  but  probably  were  too  small  to  produce  stress-strain  behavior 
conditioned  on  neighboring  grain  interactions.  This  led  to  the  focus  of 
subsequent  work  on  the  deformation  inside  individual  microplastic  grains 
and  at  their  boundaries.  Later,  careful  comparisons  of  predicted  and  mea¬ 
sured  strain  fields  outside  of  microplastic  grains  showed  that  a  common 
condition  was  for  the  matrix  surrounding  microplastic  grains  to  be  elas¬ 
tic.  The  external  elastic  strain  magnitudes  were  also  used  to  calculate 
the  depth  of  the  microplastic  strains  within  individual  grains.  The 
deeper  the  plasticity,  the  more  slowly  the  elastic  fields  fall  off  with 
distance  from  the  grain  boundary.  These  comparisons  gave  values  for  plas¬ 
ticity  deformation  depths  which  were  typical  of  the  average  grain  depth, 
suggesting  that  the  microplastic  deformation  extended  to  and  was  also 
interrupted  by  the  subsurface  boundaries. 

A  numerical  procedure  was  developed  to  calculate  from  Eshelby's  solutions 
the  correct  stress  and  strain  fields  inside  and  outside  of  a  microplastic 
grain  in  the  presence  of  a  free  surface.  Results  of  sample  calculations 
show  that  if  the  plastic  deformation  within  a  grain  is  uniform  and  volume 
conserving,  the  effect  of  the  "free  surface"  correction  is  too  small  to  be 
seen  with  our  current  experimental  sensitivity.  However,  a  singularity  in 
the  deformation  is  found  at  the  boundary  at  the  surface.  This  must  be 
relaxed  by  plastic  flow.  This  condition  may  modify  the  boundary  con¬ 
straint  and  make  itself  visible  through  the  indirect  effect  of  the  changed 
constraint.  For  example,  we  have  assumed  that  there  is  no  sliding  at  the 
boundary  between  the  microplastic  grain  and  the  external  elastic  matrix. 

If  the  singularity  promotes  sliding,  strains  near  the  boundary  would  not 
be  correctly  given  by  the  Eshelby  model. 
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An  unanticipated  result  of  strain  measurements  near  the  boundaries  of 
microplastic  grains  was  the  discovery  of  sites  in  which  motion  about  an 
external  stress-local  strain  hysteresis  loop  was  anticlockwise.  Under 
NSF^  funding,  an  explanation  for  this  phenomenon  was  sought  in  the  con¬ 
sideration  of  a  dual-domain  model  in  which  a  second  region  of  uniform 
deformation  was  embedded  at  the  boundary.  The  boundary  domain  was  allowed 
to  have  mechanical  properties  different  than  the  remainder  of  the  grain. 
Analytic  expressions  for  some  of  the  simpler  stress-strain  behaviors  of 
this  composite  structure  were  derived.  These  showed  that  with  appropriate 
parameters  a  dual  domain  structure  could  produce  the  observed  anticlock¬ 
wise  motion  at  the  boundary,  while  deformation  in  the  center  of  the  grain 
retained  the  normal  clockwise  motion.  The  boundary  behavior  is  the  result 
of  stress  fields  generated  by  constraint  of  the  boundary  domain  by  the 
rest  of  the  microplastic  grain.  Under  NAVAIR/NADC  funding  a  numerical 
procedure  was  developed  to  calculate  the  stress-strain  response  for  arbi¬ 
trary  material  parameters  and  local  yield  strengths.  By  comparing  results 
of  calculated  stress-strain  behavior  of  this  model  with  experimental  data, 
flow  stresses  were  determined  for  A1  2219-T851  as  a  function  of  location 
within  a  large  microplastic  grain.  These  values  were  much  smaller  than 
the  bulk  yield  strength  and  nearly  uniform  within  the  grain.  Near  the 
boundaries  and  inside  the  grain,  an  abrupt  increase  in  flow  stress  was 
found,  and  just  outside  the  grain,  the  strains  were  indistinguishable  from 
that  resulting  from  elastic  deformation. 

The  early  "dual -domai n"  model  was  derived  using  a  von  Mises  criterion  for 
a  uniaxial  stress  state  or  "stationary"  plastic  deformation  criterion. 

The  calculation  was  also  uniaxial  -  the  transverse  fields  created  by 
constraint  of  the  microplastic  grain  being  neglected.  This  model  lacks  a 
Bauschinger  effect  internal  to  the  material.  In  fully  reversed  loading, 
the  Bauschinger  response  predicted  stems  from  a  reaction  stress  caused  bv 
constraint  of  the  microplastic  grain.  However,  with  arbitrary  load 
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sequences,  the  absence  of  an  internal  stress  state  in  the  material  itself 
becomes  quite  apparent.  To  increase  modeling  flexibility,  a  numerical 
procedure  was  developed  to  calculate  stress-strain  response  which  included 
the  biaxial  field  and  for  a  Prager  or  "kinematic"  yielding  criterion. 

Also,  simple  analytic  expressions  for  some  of  the  deformation  behaviors  of 
the  models  were  derived  and  examined.  These  included  descriptions  for  the 
stress-strain  hysteresis  loop  widths  for  constant  amplitude  and  for  step 
changes  in  external  stress  range.  It  was  found  that  transient  variations 
in  loop  width  with  a  stress  range  change  were  intrinsic  to  the  "station¬ 
ary"  model,  but  absent  beyond  the  second  cycle  after  stress  range  change 
with  the  "kinematic"  model.  For  the  "stationary"  model  a  procedure  was 
devised  to  calculate,  from  experimental  data,  local  values  of  flow  of 
stress  (eg),  plastic  modulus  (EP),  and  the  depth  of  the  microplastic 
deformation  within  a  grain.  Application  involved  the  obtaining  of  loop 
widths  for  a  prescribed  series  of  loading  sequences,  including  determining 
the  consequence  of  change  in  stress  range  as  a  function  of  the  new  stress 
amplitude.  In  concurrence  with  the  dual  domain  results,  values  found  by 
this  technique  for  a  large  (330  pm)  grain  after  500  fatigue  cycles  were  Og 
~  225  MPa  and  EP  =  0  MPa.  This  should  be  compared  to  a  bulk  cycle  yield 
strength  for  A1  2219-T851  of  360  MPa. 

5.  Options  to  replace  the  dual  domain  model  of  deformation  of  a  microplastic 
grain  have  been  explored.  Numerical  implementation  of  a  procedure  to  cal¬ 
culate  a  physically  accurate  three-dimensional  distribution  of  strains  and 
stresses  within  a  grain  is  practical.  This  would  give  a  much  better  in¬ 
sight  into  the  mechanism  of  evolution  of  the  deformation,  especially  at 
the  grain  boundaries.  However,  an  analysis  showed  us  that  code  could  not 
be  developed  within  the  scope  of  this  program  and  this  activity  was 
terminated . 

6.  As  a  precursor  to  a  selection  of  the  correct  criterion  for  crack  initia¬ 
tion  within  a  grain,  the  behavior  of  the  hysteretic  energy  criterion  was 
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examined  in  terms  of  its  response  to  material  parameters  such  as  local 
flow  stress  and  plastic  modulus.  We  find,  for  example,  that  in  the  "sta¬ 
tionary"  deformation  model  the  energy  deposited  into  the  grain  per  cycle 
goes  down  if  the  flow  stress  drops  below  a  critical  value.  If  accumulated 
hysteretic  energy  is  driving  the  initiation,  then  the  response  of  crack 
initiation  rates  to  certain  load  sequences  will  have  characteristic  behav¬ 
iors  which  should  make  identification  of  the  energy  criterion  as  the 
correct  criterion  relatively  straightforward. 

6  .2  Conclusions  About  Def  ormat  i  on  Behavior  of  A1  2219-T851 

1.  With  the  exception  o'  ve^y  la-ge  (>  350  pm)  grains  at  high  stress  ampli¬ 
tudes  (>  85*  yield',  microplastic  deformation  is  highly  localized  and  con¬ 
tained  within  individual  surface  grains.  This  is  confirmed  by  an  equality 
between  the  elastic  strain  field  predictions  of  Eshelby’s  model  for  sites 
external  to  a  deformed  grain  and  experimental  observation. 

2.  Deformation  strains  within  a  grain  can  be  treated  as  uniform  except  near 
the  bounda-ies  where  anomalous  stress-strain  response  has  been  found.  The 
boundary  constraint  condition  which  leads  to  an  anticlockwise  motion  about 
the  external  stress-local  strain  loop  in  a  grain  near  the  boundary  is  not 
completely  understood.  However,  we  have  shown  that  a  sufficient  condition 
for  the  behavior  is  a  slightly  higher  flow  stress  near  the  boundaries 
compared  to  the  grain  interior.  A  full  understanding  of  this  behavior  is 
needed  only  for  those  materials  in  which  crack  initiation  at  the  grain 
boundaries  is  an  important  process.  However,  this  phenomena  may  influence 
the  propagation  of  small  cracks  across  grain  boundaries. 

3.  Values  of  flow  stress,  plastic  modulus  and  depth  of  deformation  have  been 
calculated  for  selected  grains.  More  needs  to  be  done  to  define  the  evo¬ 
lution  of  these  properties,  but  presently,  the  following  scenario  appears 
likely.  In  the  early  stages  of  fatigue,  only  a  very  thin  (~  2  urn)  layer 
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of  the  surface  is  microplastic.  This  depth  is  comparable  to  dislocation¬ 
cell  sizes  in  aluminum,  suggestive  of  a  connection.  The  surface  flow 
stresses  is  low  (<  150  MPa  compared  to  a  360  MPa  bulk  value),  but  the 
plastic  modulus  is  very  large  {>  10^  MPa).  The  plastic  deformation  expe¬ 
rienced  by  the  surface  is  therefore  extremely  small  and  indistinguishable 
from  elastic  except  by  the  most  sensitive  of  measurements.  With  fatigue, 
long-range  slip  develops  in  the  larger  grains  in  the  surface.  Experimen¬ 
tal  values  of  plastic  modulus  for  these  softened  grains  are  EP  ~  0.  Thus, 
the  material  in  the  grains  is  essentially  el astic-perfectly  plastic. 
Accurate  values  of  the  flow  stresses  (cQ)  for  some  grains  have  been  ob¬ 
tained  by  analysis  of  experimental  measurements  of  loop  widths  for  se¬ 
lected  load  sequences  and  load  levels.  These  give  og  values  of  approxi¬ 
mately  200  MPa  after  some  hundreds  to  thousands  of  cycles  of  fatigue  at 
stress  amplitudes  well  below  the  bulk  elastic  limit.  Less  accurate 
estimates  of  flow  stress  based  upon  a  single  loop  width  measurement  have 
been  derived  vs  grain  size,  under  the  assumption  that  Ep  =  0  for  all 
grains.  This  analysis  shows  a  flow  stress  decreasing  with  grain  size,  and 
uncovers  examples  of  grains  in  which  the  apparent  og  after  fatigue  is  less 
than  100  MPa. 

6 .3  Conclusions  About  Initiation  Lifetime  Modelling 

1.  Microplastic  deformation  which  leads  to  fatigue  crack  initiation  can  be 
highly  localized,  being  confined  to  individual  large  surface  grains. 

2.  Because  of  the  constraint  of  such  localized  plasticity  by  an  otherwise 
elastic  surface,  stresses  within  a  grain  can  depart  markedly  from  that 
externally  appl ied . 

3.  For  an  A1  2219-T851  alloy,  a  progressive  fatigue-induced  softening  of  a 
microplastic  grain  occurs.  Local  flow  stresses  apparently  can  fall  to 
less  than  30%  of  bulk  cyclic  yield  strength  values.  The  plastic  moduli  of 
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such  grains  are  small  and  the  deformation  is  nearly  elastic  -  perfectly 
pi astic . 

The  strains  within  a  grain  which  control  initiation  can  be  calculated,  at 
least  within  the  grain  interiors,  with  reasonable  reliability.  These 
respond  to  the  external  load  spectra  in  a  fashion  determined  by  the  defor 
mati on-induced  reaction  stress,  in  turn,  established  by  the  local  flow 
stress.  Herein  is  an  obvious  source  of  load  amplitude  and  sequence 
effects  on  fatigue  lifetime. 

A  complete  procedure  to  calculate  lifetimes  on  the  basis  of  such  a  defor¬ 
mation  model  awaits  the  development  of  a  model  to  relate  progressive 
changes  in  flow  stresses  within  a  grain  to  the  stress-strain  history  it 
experiences.  Candidate  criteria  to  describe  the  local  onset  of  yielding 
should  still  be  critically  reviewed  to  refine  the  accuracy  of  predicted 
stress-strain  behaviors  during  load  transients.  Finally,  the  identifica¬ 
tion  of  the  controlling  criterion  for  crack  initiation  must  be  made  from 
several  competing  candidates.  This  criterion  would  be  used  to  identify 
the  time  of  initiation  from  predicted  local  strains  in  response  to  an 
applied  spectra. 


A  capability  to  predict  the  stochastic  initiation  behavior  of  such  a 
deformation  model  would  entail  development  of  numerical  methods  to  collate 
the  response  of  m.jltiple  potential  crack  initiation  sites. 
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DEFENSE  TECHNICAL  INFORMATION  CENTER  (DTIC),  Bldg. #5,  Cameron  Station 
Alexandria,  VA  22314 

(Attn:  Administrator).  ....... 

FAA,  Washington,  D.C.  20591 

(Attn:  J.  R.  Soderquist,  AW-103).  ..... 

FAA,  Technical  Center,  Atlantic  City,  NJ  08405 

(Attn:  L.  Neri,  Code  AC~-330;  M.  Caiafa,  Code  ACT-033). 

NASA  HEADQUARTERS,  Washington,  D.  C.  20546 

(Attn:  Airframes  Branch,  FS-120).  ..... 

(Attn:  OAST/RM,  Dr.  D.  Mulville).  ..... 

NASA,  George  C.  Marshall  Space  Flight  Center,  Huntsville,  AL  35812 
(Attn:  E.  E.  Engler,  S&E-ASTN-ES ) .  ..... 

(Attn:  R.  Schwinghamer ,  S&E-ASTN-M).  .... 

NASA,  Langley  Research  Center,  Hampton,  VA  23365 

(Attn:  Dr.  J.  R.  Davidson,  MS  188E;  Dr.  J.  Starnes,  MS-190; 

Dr.  M.  Mikulus,  H.  Bohan,  and  Dr.  C.  P.  Blakenship 
MS  189M).  ........ 

NASA,  Lewis  Research  Center,  Cleveland,  OH  44135 

(Attn:  Dr.  C.  Chamis,  MS  49-6;  M.  Hershberg,  MS  49-6). 
NAVAIRSYSCOM,  Washington,  D.C.  20361 

(Attn:  AIR-00D4).  ....... 

(Attn:  AIR— 530).  ........ 

(Attn:  AIR-5302D) .  ....... 

(Attn:  AIR-5302).  ....... 

(Attn:  AIR-5302F).  ....... 

(Attn:  AIR-53032D) .  ....... 

(Attn:  AIR-931B).  ....... 

NAVPGSCHL ,  Monterey,  CA  95940 

(Attn:  Prof.  R.  Ball,  Prof.  M.  H.  Bank,  Prof.  K.  Challenger). 
NAVSEASYSCOM,  Washington,  D.C.  20360 

(Attn:  C.  Zannis,  SEA-05R25)).  ..... 

NAVSEC,  Arlington,  VA  20360 

(Attn:  NSEC-6101E).  ....... 
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NAVSHIPRANDCEN ,  Annapolis,  MD  21403 

(Attn:  H.  Edlestein,  Code  2870).  ..... 

NRL,  Washington,  D.C.  20375 

(Attn:  Dr.  I.  Wolock,  Code  6122;  Dr.  C.  I.  Chang, 

and  Dr.  R.  Badaliance).  ..... 

NSWC ,  WHITE  OAK  LABORATORY,  Silver  Spring,  MD  20910 

(Attn:  Dr.  J.  Goff,  Materials  Evaluation  Branch,  Code  R-34  . 
(Attn:  Dr.  J.  M.  Augl ) .  ...... 

ONR,  800  N.  Quincy  St.,  Arlington,  VA  22217 

(Attn:  A.  Kushner  Code  432/A;  Y.  Rajapakse,  Code  1132SM) 

ONT,  800  N.  Quincy  Street,  Arlington,  VA  22217 

(Attn:  Cdr .  D.  Brown,  (OCNR-212).  ..... 

PLASTEC ,  Piratinny  Arsenal,  Dover,  NJ  07801 
(Attn:H.  Pebly). 

(Attn:  Librarian,  Code  DRDAR-SCM-0,  Bldg.  351-N). 

ARMY  MATERIALS  TECHNOLOGY  LABORATORY 
Watertown,  MA  02171 

(Attn:  D.  Oplinger ,  SLCMT-MS).  ..... 

U.  S.  ARMY  APPLIED  TECHNOLOGY  LABORATORY,  USARTL ,  (AVRADCOM), 

Ft.  Eustis,  VA  23604 

(Attn:  J.  Waller;  T.  Hazza) .  ..... 

U.  S.  ARMY  AIR  MOBILITY  R&D  LABORATORY,  Ft.  Eustis,  VA  23604 

(Attn:  H.  Reddick).  ....... 

U.  S.  ARMY  R&T  LABORATORY  (AVRADCOM),  Ames  Research  Center, 

Moffet  Field,  CA  94035 

(Attn:  F.  Immen,  DAVDL-AS-MS  207-5).  .... 

U.  S.  NAVAL  ACADEMY,  Annapolis,  MD  21402 

(Attn:  Dr.  R.  D.  Jamison,  Mechanical  Engineering  Department). 

DAVID  TAYLOR  NAVAL  SHIP  RESEARCH  &  DEVELOPMENT  CENTER, 

Annapolis,  MD  21402 

(Attn:  E.  T.  Camponeschi,  Code  2844;  R.  Crane,  Code  2844).  . 

DAVID  TAYLOR  NAVAL  SHIP  R&D  CENTER 
Bethesda,  MD  20084 

(Attn:  A.  Macander,  Code  1720).  ..... 
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NAVAIRDEVCEN,  Warminster,  PA  18974 

(Attn:  Code  8131).  .......  3 

(Attn:  Code  09L2).  .......  2 


NADC-87046-60 


DISTRIBUTION  LIST 
NON-GOVERNMENT  ACTIVITIES 

NO.  OF 
COPIES 

ANAMET  LABORATORIES,  100  Industrial  Hyw. ,  San  Carlos,  CA  94070 

(Attn:  Dr.  R.  Arnold).  .......  1 

ALCOA  DEFENSE  SYSTEMS  CORP . ,  16761  Via  delCampo  Court, 

San  Diego,  CA  92127 

(Attn:  D.  Myers).  .......  1 

AVCO,  Specialty  Materials  Div.,  2  Industrial  Avenue,  Lowell,  MA  01851 

(Attn:  Mr.  W.  F.  Grant).  ......  1 

BATTELLE  COLUMBUS  LABORATORIES,  Metals  and  Ceramics  Information  Center 

505  King  Avenue,  Columbus,  OH  43201.  ....  1 

BEECH  AIRCRAFT  CORP.,  4130  Linden  Avenue,  Dayton,  OH  45432 

(Attn:  M.  B.  Goetz).  .......  1 

BELL  AEROSPACE  COMPANY,  Buffalo,  NY  14240 

(Attn:  F.  M.  Anthony,  Zone  1-85).  .....  1 

BELL  HELICOPTER  CO.,  Fort  Worth,  TX  76101 

(Attn:  M.  K.  Stevenson).  ......  1 

BENDIX  PRODUCTS,  Aerospace  Division,  South  Bend,  IN  46619 

(Attn:  R.  V.  Cervelli).  ......  1 

BOEING  CO.,  P.  0.  Box  3707,  Seattle,  WA  98124 

(Attn:  J.  McCarty,  J.  Quinliven,  and  Dr.  R.  June).  .  .  3 

BOEING  CO.,  Vertol  Division,  P.0.  Box  16858,  Philadelphia,  PA  19143 

(Attn:  R.  L.  Pinckney).  ......  1 

(Attn:  D.  Hart).  ........  1 

(Attn:  C.  Albrecht).  .......  1 

BOEING  CO.,  Wichita,  KS  67277-7730.  .....  1 

(Attn:  J.  Avery).  .......  1 

(Attn:  R.  Waner).  .......  1 

CABOT  CORPORATION,  Billerica  Research  Center,  Billerica,  MA  01821.  1 

DEPARTMENT  OF  TRANSPORTATION,  Kendall  Square,  Cambridge,  MA  02142.  1 

(Attn:  Dr.  Ping  Tong,  DTS  76,  TSC).  ....  1 

DREXEL  UNIVERSITY,  Philadelphia,  PA  19104 

(Attn:  Dr.  P.  C.  Chou).  ......  1 

(Attn:  Dr.  A.  S.  D.  Wang).  ......  1 

E.  I.  DuPONT  COMPANY,  Textile  Fibers  Department,  Chestnut  Run  Location 
CR701,  Wilmington,  DE  19898 

(Attn:  V.  L.  Bertarelli).  ......  1 

FAIRCHILD  REPUBLIC  CO.,  Farmingdale,  L.I.,  NY  11735 

(Attn:  Mr.  Frank  Costa).  ......] 

GEORGIA  INSTITUTE  OF  TECHNOLOGY ,  Atlanta,  GA  30332 

(Attn:  (L.  Rehfield)).  .......  1 

GENERAL  DYNAMICS /CONVAIR,  San  Diego,  CA  92138 

(Attn:  Dr.  R.  Dunbar).  .......  1 

GENERAL  DYNAMICS,  Fort  Worth  Division,  PO  Box  748,  Fort  Worth,  TX  76101 
(Attn:  J.  A.  Fant).  .......  1 

(Attn:  Composite  Structures  Eng.  Dept.).  ....  1 
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NORTHROP  AIRCRAFT  CORP . ,  One  Northrop  Avenue,  Hawthorne,  CA  90250 
(Attn:  Dr-  M.  Ratwani,  B.  Butler  and  R.  Whitehead). 

PURDUE  UNIVERSITY,  School  of  Aeronautics  and  Astronautics, 

West  Lafayette,  IN  47907 

(Attn:  Dr.  C.  T.  Sun).  ...... 

PROTOTYPE  DEVELOPMENT  ASSOCIATES,  INC.,  1560  Brookhollow  Drive 
Santa  Ana,  CA  92705 

(Attn:  E.  L.  Stanton).  ...... 

ROCKWELL  INTERNATIONAL,  Columbus,  OH  43216 

(Attn:  M.  Schweiger).  ...... 

ROCKWELL  INTERNATIONAL,  Los  Angeles,  CA  90009 

(Attn:  Dr.  Lackman).  ...... 

(Attn:  W.  O'Brien).  ...... 

ROCKWELL  INTERNATIONAL,  Tulsa,  OK  74151 

(Attn:  F.  Kaufman).  ...... 

ROHR  CORP.,  Riverside,  CA  92503 

(Attn:  Dr.  F.  Riel).  ...... 

SIKORSKY  AIRCRAFT,  Stratford,  CT  06622 

(Attn:  S.  Garbo).  ...... 

J.  P.  STEVENS  &  CO.,  INC.,  New  York,  NY  10036 

(Attn:  H.  I.  Shulock).  ...... 

TELEDYNE  RYAN  AERONAUTICAL  CO.,  San  Diego,  CA  92138 

(Attn:  R.  Long).  ....... 

UNIVERSITY  OF  DAYTON  RESEARCH  INSTITUTE,  300  College  Park  Avenue, 
Dayton,  OH  45469 

(Attn:  Dr.  J.  Gallagher).  ..... 
UNIVERSITY  OF  DELAWARE,  Mechanics  &  Aerospace  Eng.  Dept., 

Evans  Hall,  Newark,  DE  19711 

(Attn:  Dr.  R.  B.  Pipes,  Dr.  J.  R.  Vinson  and  Dr.  D. Wilkins 
UNIVERSITY  OF  OKLAHOMA,  Norman,  OK  73019 

(Attn:  Dr.  C.  W.  Bert,  School  of  AMNE). 

UNIVERSITY  OF  WYOMING,  Laramie,  WY  82071 

(Attn:  Dr.  D.  Adams).  ...... 

VILLANOVA  UNIVERSITY,  Villanova,  PA  19085 

(Attn:  Dr.  P.  V.  McLaughlin).  .... 

VIRGINIA  POLYTECHNIC  INSTITUTE,  Blacksburg,  VA  24061 

(Attn:  Dr.  K.  Reifsnider).  ..... 

WASHINGTON  UNIVERSITY,  School  of  Engineering  and  Applied  Science, 
Materials  Research  Laboratory,  Campus  Box  1087, 

St.  Louis,  MO  63130 

(Attn:  T.  Hahn).  ....... 
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NORTHROP  AIRCRAFT  CORP . ,  One  Northrop  Avenue,  Hawthorne,  CA  90250 
(Attn:  Dr.  M.  Ratwani ,  B.  Butler  and  R.  Whitehead). 

PURDUE  UNIVERSITY,  School  of  Aeronautics  and  Astronautics, 

West  Lafayette,  IN  47907 

(Attn:  Dr.  C.  T.  Sun).  ...... 

PROTOTYPE  DEVELOPMENT  ASSOCIATES,  INC.,  1560  Brookhollow  Drive 
Santa  Ana,  CA  92705 

(Attn:  E.  L.  Stanton).  ...... 

ROCKWELL  INTERNATIONAL,  Columbus,  OH  43216 

(Attn:  M.  Schweiger).  ...... 

ROCKWELL  INTERNATIONAL,  Los  Angeles,  CA  90009 

(Attn:  Dr.  Lackman).  ...... 

(Attn:  W.  O'Brien).  ...... 

ROCKWELL  INTERNATIONAL,  Tulsa,  OK  74151 

(Attn:  F.  Kaufman).  ...... 

ROHR  CORP.,  Riverside,  CA  92503 

(Attn:  Dr.  F.  Riel).  ...... 

SIKORSKY  AIRCRAFT,  Stratford,  CT  06622 

(Attn:  S.  Garbo).  ...... 

J.  P.  STEVENS  &  CO.,  INC.,  New  York,  NY  10036 

(Attn:  H.  I.  Shulock).  ...... 

TELEDYNE  RYAN  AERONAUTICAL  CO.,  San  Diego,  CA  92138 

(Attn:  R.  Long).  ....... 

UNIVERSITY  OF  DAYTON  RESEARCH  INSTITUTE,  300  College  Park  Avenue, 
Dayton,  OH  45469 

(Attn:  Dr.  J.  Gallagher).  ..... 
UNIVERSITY  OF  DELAWARE,  Mechanics  &  Aerospace  Eng.  Dept., 

Evans  Hall,  Newark,  DE  19711 

(Attn:  Dr.  R.  B.  Pipes,  Dr.  J.  R.  Vinson  and  Dr.  D. Wilkins 
UNIVERSITY  OF  OKLAHOMA,  Norman,  OK  73019 

(Attn:  Dr.  C.  W.  Bert,  School  of  AMNE). 

UNIVERSITY  OF  WYOMING,  Laramie,  WY  82071 

(Attn:  Dr.  D.  Adams).  ...... 

VILLANOVA  UNIVERSITY,  Villanova,  PA  19085 

(Attn:  Dr.  P.  V.  McLaughlin).  .... 

VIRGINIA  POLYTECHNIC  INSTITUTE,  Blacksburg,  VA  24061 

(Attn:  Dr.  K.  Reifsnider).  ..... 

WASHINGTON  UNIVERSITY,  School  of  Engineering  and  Applied  Science, 
Materials  Research  Laboratory,  Campus  Box  1087, 

St.  Louis,  MO  63130 

(Attn:  T.  Hahn).  ....... 
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